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ABSTRACT 


Investigating  the  effect  of  steam  ingestion  into  an  aircraft  jet  engine  is  necessary 
to  improve  understanding  of  stall  and  surge  in  transonic  axial  compressors.  Specifically, 
to  understand  the  “pop  stall”  phenomenon  experienced  by  naval  fighter  jet  aircraft  during 
steam  catapult  launches.  Steam  leakage  from  an  aircraft  carrier  catapult  system  can  be 
ingested  into  the  intake  and  cause  stall  or  surge  in  a  jet  engine  upon  takeoff.  It  is 
important  to  understand  the  conditions  under  which  this  occurs  as  the  Navy  prepares  for 
the  fielding  of  the  single  engine  F-35C,  the  aircraft  carrier  variant  of  the  Joint  Strike 
Fighter.  This  project  prepares  the  structure  and  instrumentation  to  investigate  the  inlet 
distortion  and  effects  of  steam  ingestion  on  a  transonic  axial  compressor.  A  compressor 
test  facility,  including  mechanical  equipment,  data  acquisition  system,  and  remote  digital 
control  system,  was  configured  to  test  a  transonic  compressor  rotor,  similar  to  what  will 
be  used  in  the  Joint  Strike  Fighter.  Rotor  inlet  and  exit  velocity  profiles  were  measured 
with  a  three-hole  probe  to  obtain  a  set  of  baseline  data  before  future  experiments. 
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I. 


INTRODUCTION 


The  F-35C  is  the  Navy  variant  of  the  Joint  Strike  Fighter  (JSF)  and  will  utilize  the 
F-135  engine,  which  uses  a  transonic  compressor.  Experiments  conducted  with  catapult 
launches  of  an  F-18,  at  Lakehurst  Naval  Air  Engineering  Station,  shown  in  Figures  1  and 
2,  demonstrated  that  current  jet  engines  are  susceptible  to  catastrophic  “pop  stall”  when 
exposed  to  excessive  amounts  of  steam.  The  Joint  Strike  Fighter  will  be  the  first  single 
engine  jet  fighter  in  recent  years  to  be  assigned  to  aircraft  carrier  based  operations.  The 
single  engine  nature  of  the  JSF  increases  the  severity  of  a  potentially  catastrophic  “pop 
stall,”  particularly  if  it  were  to  be  exposed  to  large  amounts  of  steam  leakage  during  a 
catapult  launch  aboard  an  aircraft  carrier  at  sea.  (Donelson,  2003) 


Figure  1.  Lakehurst  steam  ingestion  experiment.  F-18  approaching  steam  leak  during 

catapult  launch. 
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The  present  goal  of  the  Turbopropulsion  Laboratory  (TPL)  at  the  Naval  Post 
Graduate  School  (NPS)  is  to  evaluate  a  small-scale  transonic  compressor  during  stall  and 
surge  when  exposed  to  steam.  The  present  study  documents  the  equipment  preparation 
and  instrumentation  design  to  conduct  the  evaluation,  as  well  as  baseline  data  obtained. 


Figure  2.  Lakehurst  steam  ingestion  experiment.  F-18  “Pop  Stall”  caused  by  steam  leak 

during  catapult  launch. 

To  evaluate  the  “pop  stall”  phenomenon,  a  rotor  similar  to  that  used  in  military 
fighter  jets  was  necessary.  The  rotor  of  the  first  stage  of  a  transonic  compressor  was 
evaluated  in  the  present  study.  The  Sanger  Rotor  is  a  low  aspect  ratio  transonic 
compressor.  Nelson  Sanger  designed  the  compressor  stage  in  1996  at  the  NASA  Lewis 
Research  Center.  The  rotor  was  designed  using  CFD  techniques  and  provided  to  the 
Naval  Post  Graduate  School  for  experimental  testing.  (Sanger,  1996) 

Pressures  both  upstream  and  downstream  of  the  single  rotor  were  evaluated  by 
traversing  a  three-hole  pressure  probe  at  various  locations.  Prior  to  conducting  stall 
experiments,  a  three-hole  pressure  probe  was  calibrated  and  a  full  baseline  survey  was 


2 


completed.  Data  was  obtained  at  70,  80,  90,  and  100  percent  speed  with  mass  flow  rate 
varied  from  choke,  peak  efficiency,  and  to  stall. 

This  project  is  the  beginning  of  an  in  depth  study  of  stall  and  surge  in  transonic 
compressors.  The  current  project  prepares  the  instrumentation  and  baselines  the 
performance  of  a  transonic  rotor.  Future  projects  will  include  the  injection  of  steam 
upstream  of  the  rotor.  With  the  baseline  completed,  the  calculation  of  stall  margin  and  its 
change  as  steam  is  added  to  the  system  can  be  conducted. 
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II.  EXPERIMENTAL  FACILITY  AND  EQUIPMENT 


The  current  research  was  conducted  at  the  Turbopropulsion  Laboratory  within  the 
Department  of  Mechanical  and  Astronautical  Engineering  of  the  Naval  Postgraduate 
School. 

A.  TRANSONIC  COMPRESSOR  RIG 

The  Transonic  Compressor  Rig  (TCR)  was  originally  built  to  test  a  prototype 
transonic  compressor  stage  and  was  designed  by  Professor  Michael  Vavra.  The 
configuration  tested  is  shown  in  Figure  3  and  was  detailed  extensively  by  O'Brien  and 
Papamarkos.  (O’Brien,  2000;  Papamarkos,  2004)  The  TCR  was  set  up  to  evaluate  the 
Sanger  transonic  rotor.  The  test  compressor  was  driven  by  two  opposed  air  turbine 
stages,  which  were  supplied  with  air  by  an  Allison  Chalmers  axial  compressor,  and  was 
capable  of  testing  axial  compressors  up  to  30,000  RPM. 


Figure  3.  Transonic  compressor  rig  in  test  cell  with  inlet  piping  removed. 
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Sanger,  at  the  NASA  Lewis  Research  Center,  designed  the  compressor  rotor  as 
part  of  a  complete  fan  stage.  It  was  a  completely  Computational  Fluid  Dynamic  (CFD) 
designed  stage.  Details  of  the  design  and  design  considerations  are  contained  in  TPL 
Technical  Note  99-01.  (Sanger,  1999)  The  stage  was  manufactured  by  numerical 
machining,  specifically  for  testing  and  evaluation  at  the  NPS  Turbopropulsion 
Laboratory.  Several  experimental  studies  have  been  conducted  on  the  stage.  (Roberts, 
2003;  Gannon,  Hobson,  Shreeve,  2004)  The  rotor  profile  is  shown  on  the  TCR  in  Figure 
4. 


Figure  4.  View  of  the  rotor-only  configuration  tested. 


The  rotor  was  manufactured  from  a  high  strength  aluminum  alloy  (7075-T6)  and 
consisted  of  22  blades.  In  this  experiment,  the  rotor  was  tested  with  a  parabolic  spinner, 
which  replaced  the  conical  spinner  used  by  O’Brien.  (O’Brien,  2000;  Gannon  et  ah, 
2004) 

In  the  majority  of  previous  studies,  the  entire  stage  was  evaluated  to  verily  the 
accuracy  and  validity  of  CFD  design.  For  the  current  experiment,  the  rotor  was  evaluated 
without  the  stator  and  data  was  obtained  for  the  rotor  only,  which  was  simpler  to  test 
mechanically  and  the  future  steam  ingestion  tests  would  quantify  the  effect  on  the  rotor- 
only.  In  addition,  the  honeycomb  used  to  straighten  the  flow  downstream  of  the  stage 
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was  removed.  Design  parameters  for  the  rotor  are  reproduced  from  Sanger’s  “Design 
Methodology”  in  Table  1  below.  (Sanger,  1996) 


Table  1 .  Sanger  rotor  design  parameters. 


Parameter 

Value 

Pressure  Ratio 

1.61 

Tip  Speed 

33.02  m/s  (1300ft/s) 

Design  Speed 

27085  rpm 

Design  Mass  Flow 

7.75  kg/s  (17.051b/s) 

Specific  Mass  Flow 

170.88  kg/s-m2  (351bm/s-ft2) 

Specific  Head  Rise 

.246 

Tip  Inlet  Relative  Mach  Number 

1.28 

Aspect  Ratio 

1.2 

Hub/Tip  Radius  Ratio 

.51 

Rotor  Inlet  Ramp  Angle 

28.2 

Number  of  Rotor  Blades 

22 

Tip  Solidity  (Rotor) 

1.3 

Outside  Diameter 

0.2794m  (11  inches) 

Rotor  Diffusion  Factor  -  Tip 

.4 

Rotor  Diffusion  Factor  -  Hub 

.47 

A  compressor  map  from  “A  Transonic  Compressor  State  Part  1,”  (Gannon,  et  al., 
2004)  for  the  rotor  only  is  shown  in  Figures  5  and  6.  The  stage  design  point  is 
highlighted  for  the  100%  speed  line.  The  CFD  data  on  the  compressor  map  was  obtained 
using  SWIFT,  a  multi-block  computer  code  created  for  analysis  of  three-dimensional 
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Isentropic  Efficiency  Total -to -total  pressure  ratio 


viscous  flows  in  turbomachinery.  (Chima,  1998)  A  complete  stage  compressor  map, 
including  rotor  and  stator,  can  be  found  in  Appendix  B. 


Figure  5.  Total  pressure  ratio  compressor  map. 
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Figure  6.  Isentropic  efficiency  compressor  map. 
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A  steel  case  wall  was  used  to  mount  instrumentation  for  measurements.  The  four 
holes  labeled  in  Figure  7  were  designed  for  mounting  a  probe  actuator  for  flow 
measurements.  The  gold  band  represents  the  location  of  the  rotor,  with  holes  one  and 
two  located  upstream  of  the  rotor  and  holes  three  and  four  downstream.  Additional 
details  of  the  case  wall  can  be  found  in  Appendix  A. 


Figure  7.  Compressor  case  wall. 


Hole  locations  one  and  four  were  used  in  the  current  project.  Hole  one  was 
located  approximately  11.15cm  (4.39  inches)  upstream  of  the  leading  edge  of  the  rotor, 
while  hole  four  was  located  approximately  5.94cm  (2.34  inches)  downstream  of  the 
trailing  edge  of  the  rotor.  Hole  one  was  also  located  axially  forward  of  the  nose  cone, 
giving  access  to  a  complete  radius  of  inlet  flow  sampling,  from  center  to  case  wall. 
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In  the  TCR,  air  was  drawn  from  the  atmosphere  through  an  inlet  throttle  valve  and 
into  a  settling  chamber,  labeled  "1"  in  Figure  8.  A  46cm  diameter,  five-meter  long  inlet 
pipe  connected  the  chamber  to  the  test  compressor.  The  inlet  pipe  also  contained  a 
nozzle,  labeled  "2"  in  Figure  8,  for  flow  rate  measurement,  and  a  transition  duct  to  the 
27.94cm  diameter  case  wall.  (Papamarkos,  2004) 


B.  STEAM/INLET  DISTORTION  GENERATION 

A  Sussman  model  SVS600  steam  generator  was  installed  in  the  test  cell  adjacent 
to  the  TCR  to  provide  steam  for  the  compressor  inlet.  The  SVS600  contained  a  612kW 
boiler  and  could  provide  steam  at  a  Boiler  Horsepower  (BHP)  rating  of  62.4.  This 
corresponded  to  supply  steam  of  82.1  kgs  (1806  lbs)  per  hour  at  100C  (212F).  The 
design  pressure  of  the  boiler  is  1.034MPa  (150psig).  (Sussman,  2004)  The  steam  supply 
system  is  shown  in  Figure  9. 


Figure  9.  TCR  Steam  generation  system. 
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The  steam  generated  by  the  SVS600  will  be  injected  upstream  of  the  test 
compressor.  Steam  from  the  boiler  will  be  piped  through  a  manual  shutoff  valve  and  an 
orifice  plate  before  entering  the  TCR  test  cell.  Inside  the  test  cell,  steam  will  pass 
through  a  motor  operated  throttle  control  valve  before  mixing  with  atmospheric  air,  also 
controlled  with  a  motor  operated  throttle  control  valve.  A  fast  acting  solenoid  valve 
provides  the  capability  of  intermittently  injecting  steam  into  the  compressor  inlet.  Figure 
10  shows  the  interface  of  steam  ingestion  system  with  the  TCR  air  inlet.  The  inlet 
distortion  duct  drawing  in  Appendix  C  outlines  the  steam  injection  piping  connection  to 
the  compressor  intake  system. 


C.  FREE  JET  TEST  FACILITY 

The  free  jet  test  facility  was  used  to  calibrate  the  three-hole  probe.  A  mounting 
bracket  provided  positions  for  securing  probes  and  also  contained  a  sliding  bracket 
capable  of  pitch  movement.  The  free  jet  facility  used  air  from  the  Allison  Chalmers 
compressor  to  create  a  steady  flow  through  a  10.8cm  (4.25inch)  nozzle.  For  Mach 
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numbers  above  0.3,  a  dump  valve  on  the  compressor  was  throttled  to  adjust  the  flow.  For 
Mach  numbers  below  0.3,  a  manual  shutoff  valve  located  in  the  free  jet  cell  was 
activated.  The  free  jet  is  shown  with  the  probe  actuator  mounted  on  the  mounting  bracket 
in  Figure  11. 


Figure  1 1 .  Three-hole  probe  within  the  probe  actuator  in  the  free  jet  assembly. 
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III.  INSTRUMENTATION 


A.  PROBES 

The  primary  probe  used  in  this  experiment  was  the  three-hole  cobra  probe  shown 
below  in  Figure  12.  In  addition,  two  additional  probes  were  used  in  both  the  calibration 
and  experimentation  phases  of  this  project.  A  kiel  probe  was  used  to  detect  the 
stagnation  pressure  and  a  pitot-static  probe  was  used  to  detect  the  static  pressure  within 
the  free  jet.  The  probes  positioned  for  calibration  in  the  free  jet  test  facility  are  show  in 
Figure  13. 


Figure  12.  Three-hole  cobra  probe. 


B.  L.C.  SMITH  PROBE  ACTUATOR 

The  L.C.  Smith  probe  actuator,  shown  in  Figure  11,  was  used  to  traverse  and 
rotate  the  probe  for  calibration  in  the  free  jet  and  testing  in  the  TCR.  The  probe  actuator 
was  initially  set  up  in  the  same  manner  as  described  by  Greco;  however,  the  actuator  was 
configured  to  an  automatic  PC  controlled  system.  (Greco,  1995) 

C.  HP VEE/S C ANIV ALV E  SYSTEM 

The  first  calibration  process  used  the  Hewlett  Packard  Visual  Engineering 
Environment  (HPVEE)  and  Scanivalve  system  to  obtain  pressures  in  the  cobra  probe. 
The  process  used  a  setup  identical  to  that  outlined  by  Grossman,  and  screenshots  of  the 
post  processing  calculation  program  used  in  the  current  project  are  shown  in  Appendix  D. 
(Grossman,  1997) 
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Figure  13.  Probes  positioned  on  the  mounting  bracket  for  calibration  procedure  on  the  free- 

jet  test  assembly. 


D.  OMEGA  PX-138  PRESSURE  TRANSDUCERS 

The  Omega  PX-138  transducers  used  in  the  latter  part  of  the  experiment  were 
miniature  pressure  transducers,  one  of  which  is  shown  in  Figure  14.  The  pressure 
transducer  used  silicon  pressure  sensors,  which  were  temperature-compensated.  The 
PX138-015D5V  model  was  capable  of  a  differential  pressure  range  of  -15PSI  to  +15PSI 
with  a  linearity  and  hysteresis  error  of  ±0.1%  typical,  ±0.5%  max.  The  PX-138  was 
given  a  regulated  8VDC  power  supply  and  provided  an  output  of  2.5VDC  to  4.5VDC  for 
the  range  of  pressures  applied  in  the  experiment.  (Omega,  2005) 
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Figure  14.  PX- 138  pressure  transducer 

The  pressure  transducers  remained  calibrated  throughout  the  experiment.  Their 
pressure  to  voltage  ratio  was  nearly  linear  and  all  five  transducers  would  generally  read 
atmospheric  pressure  within  10  Pascals,  a  less  than  0.1%  error.  Flowever,  the  PX-138 
transducer  had  very  small  pinouts  and  came  with  a  CX-136-4  plastic  mating  connector. 
Maintaining  electrical  connectivity  was  occasionally  a  problem  due  to  the  fragile  nature 
of  the  pinouts  and  connectors. 

E.  PMD  1608FS  DATA  ACQUISITION  DEVICE 

The  Personal  Measurement  Device™  (PMD)  1608FS  Data  Acquisition  Device 
was  a  Universal  Serial  Bus  (USB)  2.0  (full-speed)  compatible,  analog  to  digital  (A/D) 
converter.  The  PMD  1608FS  had  a  dedicated  A/D  converter  for  each  channel  and  was 
capable  of  the  simultaneous  sampling  of  up  to  eight  channels.  Each  channel  was  capable 
of  sampling  at  50  kS/s  with  a  maximum  of  200  kS/s  for  all  channels.  The  PMD  1608FS 
was  powered  through  the  5 -volt  USB  port  and  it  was  controlled  with  a  PC  through  the 
USB  port  with  Matlab  programming.  A  detailed  block  diagram  from  the  user's  manual 
is  shown  in  Figure  15. 
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Figure  15.  PMD  1608FS  A/D  Module  block  diagram,  (from  Measurement  Computing 

“PMD  1608FS”,  2005) 


F.  USB  ERB-24  REMOTE  RELAY  CONTROLLER 

A  USB  ERB-24  Remote  Relay  Controller  was  used  to  control  the  traverse  of  the 
L.C.  Smith  gear.  It  was  capable  of  being  addressed  by  the  PC  via  the  USB  2.0  port.  The 
USB  ERB-24  was  an  electromechanical  relay  module  with  24  electromechanical  single 
pole  double  throw  relays.  The  module  was  configured  in  two  banks  of  eight  and  two 
hanks  of  four  relays.  Each  relay  was  capable  of  switching  6  amps  at  28VDC  with  a 
minimum  close  time  of  10  milliseconds  and  a  minimum  open  time  of  5  milliseconds. 
The  module  used  a  9VDC  power  supply  and  was  also  controlled  through  the  USB  port 
with  Matlab  programming.  A  detailed  block  diagram  from  the  user's  manual  is  shown  in 
Figure  16. 
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Figure  16.  USB  ERB24  Relay  Module  block  diagram,  (from  Measurement  Computing 

“USB-ERB24”,  2005) 


G.  MATLAB  DRIVER 

Matlab  Release  13,  Version  6.5  was  used  for  programming  and  control  of  data 
acquisition  devices.  Matlab  was  able  to  address  the  relays  and  guide  the  probe  actuator 
through  feedback  controlled  rotations.  In  addition,  Matlab  was  able  to  collect  and 
process  pressure  data  from  the  transducers  prior  to  exporting  it  to  an  Excel  spreadsheet. 
A  single  Graphical  User  Interface  (GUI)  was  created  to  perform  the  required  functions. 
The  GUI,  shown  in  Figure  17,  was  capable  of  probe  actuator  control  and  data  collection 
in  the  calibration  procedure. 
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Figure  17.  Matlab  calibration  graphical  user  interface. 

The  calibration  control  GUI  required  the  manual  input  of  the  stagnation 
temperature  of  the  jet,  read  from  the  HPVEE  system,  and  the  starting  angle,  ending  angle, 
and  number  of  calibration  data  points,  or  angles,  as  input.  The  function 
"calibrationrun.m"  was  called  from  the  Excel  control  box  and  outputted  the  required 
angles  to  use  for  calibration  and  used  the  feedback  and  control  function  "rotate. m",  which 
accessed  the  both  the  USB  relay  and  USB  data  acquisition  (DAQ)  device  to  position  the 
probe  actuator  for  each  angle. 

At  each  calibration  angle,  the  probe  would  stop  and  the  "outputtoexcel.m" 
function  called  the  "pressurecalibration.m"  function  to  take  pressure  readings  and 
calibrate  them.  The  pressure  data  was  combined  with  the  stagnation  temperature  of  the 
jet,  from  the  GUI  input,  to  calculate  the  calibration  constants.  The  raw  data  along  with 
the  processed  data  was  exported  to  an  Excel  spreadsheet.  A  detailed  block  diagram 
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showing  both  probe  the  actuator  control  and  the  data  acquisition  process  for  the  probe 
calibration  procedure  for  Matlab  programming  is  shown  in  Figure  18. 
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Figure  18.  Matlab  calibration  code  block  diagram 


The  Matlab  data  acquisition  GUI  was  adapted  from  the  calibration  GUI.  The  GUI 
was  changed  to  automatically  control  both  the  angle  and  the  traverse  of  the  probe 
actuator.  The  Excel  control  box  was  changed  to  require  input  for  the  starting  point  and 
ending  point  for  probe  traversing,  along  with  the  number  of  points  required  for  sampling. 
Once  the  data  was  entered,  the  "SAMPLE  DATA"  pushbutton  sent  the  probe  into  an 
automatic  cycle  of  probe  traverse,  rotation,  and  data  collection.  A  screenshot  of  the 
Matlab  data  collection  GUI  is  shown  in  Figure  19. 
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Figure  19.  Matlab  data  collection  graphical  user  interface. 


On  the  data  acquisition  GUI,  when  the  "SAMPLE  DATA"  pushbutton  was 
depressed,  the  "sampledatarun.m"  function  processed  the  information  that  was  inputted 
in  the  GUI  and  controlled  a  series  of  traverses  to  each  radial  position  for  data  collection. 
At  each  radial  position,  the  probe  would  sample  the  flow  field  and  process  the  pressure 
data  using  the  constants  previously  calculated  from  the  pressure  calibration  to  determine 
the  flow  relative  to  the  pressure  probe.  When  necessary,  the  "sample  data  run.m" 
function  would  input  a  new  angle  to  the  "rotate. m"  function  to  perform  a  null  yaw 
routine.  The  data  would  again  be  sampled  and  processed.  The  raw  pressure  and  probe 
actuator  angle  data,  along  with  processed  data  including  the  Mach  number  sensed  and 
flow  angle,  were  outputted  to  an  Excel  Spreadsheet.  The  Matlab  experimental  DAQ  and 
feedback  control  block  diagram  is  shown  in  Figure  20. 
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Figure  20.  Matlab  experimental  data  acquisition  and  feedback  control  code  block  diagram. 
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IV.  CALIBRATION 


A.  PURPOSE 

The  purpose  of  calibrating  a  multi-hole  pressure  probe  was  so  that  it  could  be 
placed  in  an  unknown  flow  field  to  determine  the  magnitude  and  direction  of  the  flow.  In 
this  experiment  a  three-hole  pressure  probe  was  calibrated  by  inserting  it  into  a  known 
flow  field  and  rotating  it  through  a  spectrum  of  specific  yaw  and  pitch  angles.  The  data 
from  each  position  angle  and  velocity  was  converted  into  a  calibration  database,  utilizing 
a  set  of  non-dimensional  coefficients.  Accurate  predictions  of  the  velocity  vector  in  an 
unknown  flow  field  could  then  be  determined  by  interpolating  the  non-dimensional 
coefficients.  (Johansen,  Rediniotis,  Jones,  2001) 

B.  OVERVIEW 

The  calibration  experiment  was  conducted  twice.  In  each  case,  the  three-hole 
cobra  probe  was  calibrated  using  the  free  jet  test  rig.  A  pitot-static  probe,  a  stagnation 
probe,  and  a  temperature  probe  were  used  to  collect  the  calibration  data.  The  probes 
were  mounted  on  a  bracket  and  situated  around  the  free  jet.  The  three-hole  probe  was 
mounted  in  the  L.C.  Smith  actuator,  capable  of  both  axial  and  rotational  movement. 

Pressure  data  was  initially  acquired  using  the  Scanivalve  system.  The  Scanivalve 
system  required  an  excessive  amount  of  time,  as  it  contained  only  one  pressure  transducer 
and  had  to  step  through  each  pressure  and  then  wait  for  tube  pressure  to  stabilize  before 
taking  a  reading. 

Once  an  initial  calibration  run  was  completed,  the  setup  was  converted  to  a 
computer  logic  controlled  system  using  USB  devices  accessed  with  Matlab.  Data  was 
acquired  through  these  devices  and  the  run  was  completed.  The  response  time  of  the 
computer  feedback  control  system  greatly  reduced  the  required  time  for  the  three-hole 
probe  calibration. 

In  addition,  the  Scanivalve  system  was  replaced  with  five  analog  output  pressure 
transducers.  Each  pressure  transducer  was  connected  to  the  USB  data  acquisition  device. 
The  data  from  the  pressure  transducers  was  verified  to  be  accurate  when  compared  to  the 
previous  data  obtained  with  the  Scanivalve.  In  both  cases  the  resulting  calibration 
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coefficients  for  the  pitch  angles  were  not  sufficiently  discriminatory  and  pitch 
measurements  were  found  to  be  unreliable,  hence  pitch  was  not  used  in  the  experimental 
data  acquisition  for  the  three-hole  probe. 

C.  CALIBRATION  PROCEDURE 

The  probe  actuator  containing  the  three-hole  probe  was  attached  to  the  free  jet  test 
assembly  previously  shown  in  Figure  12.  The  jet  nozzle  had  a  diameter  of  10.8cm 
(4.25inches)  and  the  pitot-static  and  stagnation  probes  were  both  placed  at  the  same  radial 
distance  from  the  center  of  the  flow  2.86cm  (1.125inches). 

The  three-hole  probe  was  placed  into  a  Mach  0.5  flow  and  rotated  until  a 
manometer  showed  equal  pressures  on  ports  2  and  3  shown  in  Figure  2 1 .  The  center  hole 
was  designated  "1"  in  accordance  with  normal  industry  standards.  When  facing  into  the 
flow,  the  left  side  was  designated  "2"  and  the  right  side  was  designated  "3." 


Figure  2 1 .  Three-hole  probe  numbered. 

Once  the  manometer  showed  equal  pressures,  the  pressure  tubes  were  removed 
from  the  manometer  and  connected  to  the  Scanivalve  ports  5  and  7.  The  Scanivalve 
system  verified  ports  2  and  3  read  equal  pressures  on  the  HPVEE  system.  The  angle  at 
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which  the  outer  pressures  were  equal  was  designated  the  "zero"  angle  and  used  as  a 
reference  point  for  calibration. 

A  post-processing  program  was  created  using  the  HPVEE  system  to  detennine  the 
calibration  constants  for  the  three-hole  probe.  Calibration  data  was  obtained  using  the 
HPVEE/Scanivalve  system  described  by  Grossman.  (Grossman,  1997)  Screen  captures 
of  the  HPVEE  system  are  shown  in  Appendix  D. 

In  the  free  jet  test  assembly,  the  pitot-static  and  stagnation  probes  were  used  to 
calculate  the  Mach  number.  The  Mach  number  was  increased  and  decreased  by  adjusting 
the  mass  flow  rate  through  the  nozzle.  This  was  conducted  by  throttling  the  dump  valve 
on  the  Allison  Chalmers  compressor  for  Mach  numbers  above  0.3  and  by  throttling  the 
control  valve  in  the  cell  for  below  Mach  0.3. 

The  calibration  was  conducted  by  rotating  the  probe  from  -30°  to  30°  with  a 
manual  electric  control  switch.  Pressure  and  temperature  readings  were  taken  at  five- 
degree  increments  and  exported  to  an  Excel  database.  Each  complete  rotation  was 
conducted  for  Mach  numbers  0.85,  0.7,  0.5,  0.3,  .1  and  0.05  with  the  probe  perpendicular 
to  the  flow. 

D.  CALLIBRATION  PROCEDURE  WITH  AUTOMATED  CONTROLS  AND 

INDIVIDUAL  PRESSURE  TRANSDUCERS 

The  probe  actuator  was  adapted  to  a  feedback  control  system  using  the  USB 
controlled  relays,  and  a  USB  A/D  converter.  The  feedback  was  provided  through  the 
Measurement  Computing  (MC)  Analog  Input  device  PMD-1608FS,  which  measured  the 
voltage  output  of  the  potentiometers  on  the  probe  actuator  and  was  accessed  by  Matlab 
via  a  USB  cable.  Control  of  the  actuator  was  also  provided  via  USB  cable,  with  Matlab 
accessing  relays  on  the  MC  ERB-24  to  rotate  and  traverse  the  probe.  Matlab  control 
programming  can  be  found  in  Appendix  G.  The  feedback  control  system  and  data- 
acquisition  wiring  diagram  is  shown  in  Figure  22. 
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Figure  22.  Probe  actuator  and  data  acquisition  wiring  diagram. 
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Omega  Engineering  PX-138  Pressure  transducers  were  integrated  into  the 
pressure  lines  to  replace  the  Scanivalve  system.  Five  pressure  transducers  were  used  to 
provide  individual  pressure  outputs  for  the  three-hole,  static,  and  stagnation  pressure 
probes.  Each  pressure  transducer  produced  an  analog  signal  (voltage),  which  was 
digitized  by  the  PMD-1608FS  and  converted  to  engineering  units  (Pascals)  using  a 
Matlab  calibration  program.  The  Matlab  program  also  calculated  the  probe  calibration 
constants  and  exported  all  of  the  data  to  an  Excel  Spreadsheet. 

The  transducers  were  calibrated  by  creating  a  table  of  voltage  outputs  for  known 
pressures.  The  transducers  had  two  ports  to  measure  a  differential  pressure.  One  port  of 
each  transducer  was  connected  to  a  pressure  standard.  To  verify  accuracy,  the  calibration 
pressure  was  then  connected  to  the  other  pressure  port  on  each  transducer  and  the  outputs 
were  once  again  tabulated.  Once  complete,  the  pressure  transducers  were  disconnected 
from  the  calibration  pressure  and  connected  to  the  pressure  probes.  The  response  of  the 
transducer  was  found  to  be  nearly  linear  for  both  sides,  positive  and  negative  pressures.  A 
linear  interpolation  method  in  Matlab  was  used  to  find  the  calibration  for  pressures. 
When  the  PX-138  was  used  for  calibration,  it  was  commanded  to  take  30  samples  at 
500hz,  average  the  results,  and  repeat  this  process  four  times.  The  four  results  were  again 
averaged  to  create  the  output  of  each  pressure  transducer.  Calibration  data  and  Matlab 
programming  for  the  PX-138  transducers  can  be  found  in  Appendix  F. 

The  calibration  experiment  was  repeated  using  the  control  system  and  PX-138 
pressure  transducers.  Using  a  feedback  loop  for  angle  position  and  automatically  rotating 
through  the  -30°  to  30°  the  time  required  for  the  experiment  was  greatly  reduced.  Data 
obtained  for  calibration  was  comparable  to  that  obtained  by  the  Scanivalve  system. 
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E.  CALIBRATION  EQUATIONS 

1.  Purpose 

The  purpose  of  the  calibration  calculation  was  to  determine  the  non-dimensional 
variables  for  each  angle.  The  total  pressure  coefficient  (At_probe),  and  the  static  pressure 
coefficient  (As_probe)  had  to  be  determined  for  each  angle  in  a  flow  field  of  a  known 
velocity  to  be  later  used  to  determine  the  flow  angle  in  an  unknown  field.  A  sample 
calculation  of  the  calibration  constants  is  shown  in  Appendix  E  to  demonstrate  the 
calibration  process.  The  calibration  calculations  were  adapted  from  five-hole  probe 
calibration  calculations.  (Johansen,  et  al.,  2001) 

2.  Setup  and  Nomenclature 

The  center  hole  was  designated  "1"  in  accordance  with  standard  industry 
designation.  When  facing  into  the  flow,  the  left  side  was  designated  "2"  and  the  right 
side  was  designated  "3."  Pitch  and  yaw  angle  definitions  are  shown  in  Figure  23. 
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Table  2.  Calibration  nomenclature. 


Symbol 

Definition 

P 

Yaw  Angle 

Pi 

Center  Probe  Pressure 

P2 

Port  or  Left  Side  Probe  Pressure  when  facing  flow 

P3 

Starboard  or  Right  Side  Probe  Pressure  when  facing  flow 

Poj 

Jet  Stagnation  Pressure 

Psj 

Jet  Static  Pressure 

TP 

Jet  Stagnation  Temperature 

Y 

Specific  Heat  Ratio 

cP 

Constant  Pressure  Specific  Heat 

q_probe 

Probe  Pseudo  Dynamic  Pressure 

VJet 

Jet  Velocity 

Cp_i 

Left  High  Angle  Flow  Discriminator 

Cp_2 

Center  Low  Angle  Flow  Discriminator 

Cp_3 

Right  High  Angle  Flow  Discriminator 

At  Probe 

Probe  Total  Pressure  Coefficient 

As  Probe 

Probe  Static  Pressure  Coefficient 

MachNo 

Mach  Number  Calculated  with  Jet  Static  and  Stagnation  Pressure 

Probes 
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3.  Calibration  Equations 

The  pseudo  dynamic  pressure,  q_probe,  provided  a  single  pressure  value  for  each 
data  point  for  the  probe  and  was  given  as: 


q_probe  :=  P  ^ 


(p2  +  r3) 

2 


4.1 


The  jet  velocity  was  calculated  with  the  static  and  stagnation  probes  and  is  given 
as: 


VJet  := 


(y-i) 

f  p  A  Y 

Sj  1 

1  - 

•  /  2-C_-T_ 

vPojJ 

'■c 

4.2 


The  angle  flow  discriminators,  Cpj,  Cp  2,  and  Cp  3,  provided  a  sense  of  the 
magnitude  of  the  direction  of  the  flow,  relative  to  the  probe. 

r  .  (pi-pd 

p-  q_probe  43 

r  .  (Pl-P3) 

p-  q_probe  44 


(p  2  - p  3) 

q_probe 


4.5 
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The  total  and  static  pressure  coefficients  provided  the  measurable  data  to  be  used 
during  flow  field  surveys. 


(P  ■  -  P  Oj) 

LProbe  • 


"l  s  Probe  ' 


q_probe 

(P  oj  “  P  sj) 


4.7 


The  Mach  number  was  also  calculated  from  the  static  and  stagnation  probe 
readings: 


MachNo  := 


(y-i) 


(  P 


Y 


-  1 


ro 

vY  -  U 


4.8 


The  calibration  constants  for  each  speed  and  angle  were  tabulated  in  an  Excel 
spreadsheet.  When  the  table  was  complete,  the  data  was  stored  and  placed  in  a  binary 
Matlab  format  (MAT)  file  for  use  in  subsequent  flow  field  surveys. 

E.  PREDICTION  OF  VELOCITY  AND  YAW  ANGLE 

To  determine  the  unknown  velocities  and  yaw  angles,  the  Matlab  function 
GRIDDATA3  was  used.  The  calibration  data  was  accessed  through  a  saved  Matlab 
workspace  text  MAT  file  and  loaded  each  time  the  flow  field  survey  was  conducted. 
GRIDDATA3  interpolated  hyper-surfaces  using  tessellation-based  linear  interpolation 
based  on  a  Delaunay  triangulation  of  the  data.  (The  Mathworks  Inc,  2002)  Interpolation 
constants  and  Matlab  programming  are  contained  in  Appendix  I. 
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F.  COMPARISON  OF  SCANIVALVE  AND  PX-138  DATA 

The  Scanivalve  and  PX-138  systems  created  similar  calibration  databases. 
However,  the  PX-138  system  was  much  quicker,  and  produced  smoother  calibration 
curves  for  the  As  and  At  calibration  constants,  especially  in  the  lower  Mach  numbers. 
The  Scanivalve  system  produced  data  with  more  variations  when  the  flow  was  below 
Mach  0.5.  Three-dimensional  graphs  of  the  Scanivalve  and  PX-138  calibration  surfaces 
are  shown  in  Figures  24-27. 

The  calibration  curves  for  the  transducer  data  are  visually  smoother  than  that  of 
the  Scanivalve  system.  However,  in  comparing  the  calibration  surface  graphs, 
consideration  must  be  given  to  the  fact  that  the  Scanivalve  system  was  calibrated  to  a 
Mach  number  of  0.05  and  the  PX-138  Transducers  were  only  calibrated  to  a  Mach 
number  of  0.066.  The  data  was  not  calculated  as  low  for  the  PX-138  transducers  because 
the  Scanivalve  data  was  inconsistent  below  Mach  0.1  and  it  was  not  believed  the 
transducers  could  be  calibrated  for  very  low  flow  fields  as  well.  Further  calibration 
details,  including  calibration  calculations  and  variables,  software  programs  used, 
calibration  data,  and  calibration  curves,  are  shown  in  Appendices  E-I. 
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As 


Figure  24.  Scanivalve  "As"  Calibration  Surface. 
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Figure  25.  PX-138  Transducer  "As"  Calibration  Surface. 
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Figure  26.  Scanivalve  "At"  Calibration  Surface. 


Figure  27.  PX-138  Transducer  "At"  Calibration  Surface. 
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V.  TEST  RUN  PROCEDURE 


A.  OVERVIEW 

The  L.C.  Smith  probe  actuator  containing  the  three-hole  probe  was  attached  to  the 
steel  case  wall  of  the  Transonic  Compressor.  Test  runs  were  conducted  on  both  upstream 
and  downstream  positions.  The  probe  was  inserted  into  the  unknown  flow  at  various 
depths  and  data  was  collected  at  each  position.  The  process  was  repeated  for  different 
speeds  and  mass  flow  rates. 

The  data  collected  is  the  baseline  data  to  determine  the  normal  operational 
parameters  of  the  rotor  at  various  speeds  and  mass  flow  rates.  From  this  data, 
experiments  following  can  calculate  the  effects  of  steam  ingestion  on  stall  margins  and 
rotor  perfonnance. 

B.  TYPICAL  TEST  PROCEDURE 

The  upstream  position  utilized  hole  one  while  the  downstream  position  used  hole 
four.  The  L.C.  Smith  actuator  was  firmly  attached  to  the  mounting  bracket  on  the  steel 
case  wall  and  the  probe  was  controlled  through  Matlab  programming  with  a  remote  PC  as 
described  in  Chapter  III. 

For  the  upstream  position,  15  data  points  were  taken.  The  Matlab  program  used  a 
cosine(7t/2)  function  to  identify  the  points  to  sample  data.  The  cosine(7t/2)  function  was 
projected  on  to  the  distance  between  center  and  case  wall  so  that  the  frequency  of  data 
points  increased  as  the  probe  approached  the  case  wall.  The  probe  would  begin  at  the 
center  of  the  flow  cylinder  and  incrementally  traverse  towards  the  case  wall.  The 
increased  number  of  data  points  closer  to  the  case  wall  improved  the  resolution  of  the 
boundary  layer.  Positions  of  the  cosine(7i/2)  function  are  shown  in  Table  3. 

For  the  downstream  position,  15  data  points  were  also  taken.  For  this  position,  a 
cosine(7t)  function  was  used  to  sample  data  more  frequently  closer  to  the  hub  and  case 
wall.  A  depth  micrometer  was  used  to  identify  the  exact  distance  from  the  hub  to  the 
case  wall.  The  probe  would  begin  at  an  approximate  1.5mm  offset  from  the  hub,  and 
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traverse  towards  the  case  wall.  The  cosine(7t)  function  was  projected  on  to  the  distance 
between  hub  and  case  wall.  Positions  of  the  cosine(Ti)  function  are  also  shown  in  Table  3 
below. 


Table  3.  Probe  sampling  positions 


Position 

Upstream  Cos(  n/2)  in  cm 

Downstream  Cos(ti)  in  cm 

1 

13.97 

4.5 

2 

12.41 

4.44 

3 

10.86 

4.28 

4 

9.36 

4.01 

5 

7.91 

3.65 

6 

6.54 

3.23 

7 

5.26 

2.75 

8 

4.09 

2.25 

9 

3.05 

1.75 

10 

2.14 

1.27 

11 

1.38 

.85 

12 

.78 

.49 

13 

.35 

.22 

14 

.09 

.06 

15 

0 

0 

Null  yawing  was  used  in  the  downstream  position.  The  angle  of  flow  varied  from 
the  hub  to  the  case  wall  and  with  speed  and  mass  flow  rate.  To  ensure  the  accuracy  of  the 
measurement,  an  additional  loop  in  the  Matlab  programming  was  created.  Logic  was 
added  so  that  if  the  measured  angle  of  the  flow  were  greater  than  5°  from  the  angle  of  the 
probe,  the  probe  actuator  would  turn  the  probe  into  the  flow  to  null  yaw  the  probe.  The 
calibration  of  the  probe  was  assumed  to  be  very  accurate  within  the  near  linear  area 
between  ±15°  and  satisfactory  results  were  obtained.  However,  the  final  data  point 
closest  to  the  case  wall  occasionally  read  a  large  angle  close  to  +  or  -  90°.  This  caused 
the  Matlab  program  to  fail  to  output  the  final  position  as  the  angle  was  out  of  normal 
parameters.  It  was  believed  this  erroneous  data  is  caused  by  the  probe  tip's  proximity  to 
the  probe's  insertion  hole  in  the  case  wall  and  possible  swirling  of  flow  in  that  region. 

Each  test  run  was  conducted  at  a  constant  speed,  with  mass  flow  rate  adjusted  for 
each  set  of  data  points.  Runs  were  conducted  at  70,  80,  90,  and  100  percent  speed.  100 
percent  speed  corresponded  to  the  design  speed  of  27,085  revolutions  per  minute. 
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On  each  test  run,  the  first  set  of  data  was  conducted  with  the  compressor  throttle 
fully  open.  This  setting  was  labeled  as  the  "Choke"  setting.  For  the  second  run,  the  mass 
flow  rate  was  reduced  to  the  predetennined  "Peak  Efficiency"  setting  for  the  each  speed 
in  accordance  with  the  Mass  flow-Isentropic  efficiency  graph  contained  in  Appendix  B. 
The  final  run  of  each  series  was  with  a  reduced  flow  corresponding  to  the  "near  stall" 
mass  flow  rate.  Table  four  identifies  the  rotor  speeds  and  mass  flow  rates  used  for 
calculations. 


Table  4.  Rotor  speed  and  mass  flow  rates. 


Speed 

Mass  Flow  Rates  (Kg/sec) 

Speed 

Speed 

(RPM) 

Choke 

Peak  Efficiency 

Near  Stall 

70% 

18959.5 

5.95 

5.58 

4.61 

80% 

21668 

6.73 

6.47 

5.49 

90% 

24376.5 

7.39 

7.24 

6.21 

100% 

27085 

8.09 

7.97 

7.69 
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VI.  RESULTS  AND  DISCUSSION 


A.  ANALYSIS 

1.  Upstream  Data 

Inlet  surveys  were  conducted  at  70,  80,  90,  and  100  percent  speed.  The  radial 
Mach  number  distributions  are  shown  in  Figures  28-31  with  yaw  angle  and  stagnation 
pressure  in  Appendices  J  and  K. 
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Figure  30.  Inlet  Mach  number  at  90%  speed 
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2.  Downstream  Data 

Exit  surveys  were  conducted  downstream  of  the  rotor.  Results  are  shown  in 
Figures  32-39  for  rotor  speeds  of  70,  80,  90,  and  100%  with  Mach  number  and  flow 
angle  plotted  against  the  hub-to-tip  radius  ratio.  Graphs  of  the  downstream  stagnation 
pressure  ratio  for  each  speed  are  shown  in  Appendix  L. 
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Figure  32.  Exit  Mach  number,  70%  speed. 


Figure  33.  Exit  Mach  number,  80%  speed. 
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Figure  34.  Exit  Mach  number,  90%  speed. 


Figure  35.  Exit  Mach  number,  100%  speed. 
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Figure  36.  Exit  flow  angle,  70%  speed. 


0.65  0.7  0.75  0.8  0.85  0.9  0.95  1 


Hub-to-tip  radius  ratio 


Figure  37.  Exit  flow  angle,  80%  speed. 
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Figure  38.  Exit  flow  angle,  90%  speed. 


Figure  39.  Exit  flow  angle,  100%  speed. 
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3.  Diffusion  Factor 

The  diffusion  factor  is  a  method  of  assessing  the  blade  loading  in  an  axial 
compressor.  For  a  simple,  two-dimensional  geometry,  the  diffusion  factor  was  calculated 
using  the  equation: 

V  2  V02 

DF  :=  1 - + - 

V  i  2-ct  V  j 

Where  Vi  and  V2  were  the  inlet  and  exit  velocities  relative  to  the  blade,  V02  was 
the  turning  of  the  flow,  assuming  the  flow  is  axial  at  the  blade  entrance,  and  a  is  the 
solidity  of  the  blade  corresponding  to  the  flow.  It  was  believed  that  diffusion  factor 
values  above  0.6  could  indicate  blade  stall  and  values  of  0.45  may  be  a  typical  design 
choice.  (Cumpsty,  1989) 

The  relative  inlet  velocity,  Vi  was  calculated  using  the  velocity  of  the  rotor  and  by 
taking  the  area  average  of  the  velocities  found  upstream  and  calculating  the  flow  at  the 
leading  edge  of  the  rotor  assuming  the  decreasing  area  acts  as  a  nozzle.  The  relative  exit 
velocity,  V2,  was  calculated  using  velocity  of  the  rotor  from  the  downstream  probe  angle 
and  velocity  measurement  and  calculating  the  radial  and  axial  components  of  flow.  From 
these  components,  the  relative  exit  velocity  was  determined. 

The  solidity  was  calculated  using  data  from  the  rotor  blade  that  had  been  used  in 
the  SWIFT  calculations.  A  graph  of  the  solidity  can  be  found  in  Appendix  O. 

The  diffusion  factor  was  calculated  for  each  set  of  speed  data.  A  sample  set  of 
calculations  for  100%  Speed,  Peak  Efficiency,  near  the  midpoint  of  the  blade  is  shown  in 
Appendix  O.  Three-dimensional  representations  of  the  diffusion  factor  are  shown  for 
Choke,  Peak  Efficiency  and  Stall  in  Figures  40-42. 
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Figure  40.  Diffusion  Factor  at  Choke. 


Figure  4 1 .  Diffusion  Factor  at  Peak  Efficiency. 
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Figure  42.  Diffusion  Factor  at  Near  Stall. 

B.  OBSERVATION 

1.  Tip  Stall 

The  experimental  data  predicted  a  tip  stall  as  the  mass  flow  rate  is  decreased.  In 
Figures  40-42,  the  diffusion  factor  peaks  near  the  tip.  Some  research  has  been  conducted 
on  the  shock  structure  and  location  of  the  tip  leakage  vortex  originating  at  the  location  of 
the  maximum  pressure  differential  at  the  tip  of  the  blade.  (Gannon  et  ah,  2004) 

2.  Magnitude  of  Diffusion  Factor 

As  the  speed  is  increased  or  the  mass  flow  rate  is  decreased,  the  diffusion  factor 
increased.  This  demonstrated  the  increased  blade  loading  for  those  conditions.  The 
general  value  for  the  diffusion  factor  of  an  axial  compressor  is  0.45  and  the  value  for  the 
diffusion  factor  at  stall  is  0.6.  (Cumpsty,  1989)  In  Figure  42,  which  shows  the  diffusion 
factor  at  near  stall  for  speeds  70-100%,  the  diffusion  factor  approached  0.6  at  the  tip  for 
all  speeds  measured.  The  high  diffusion  factors  at  the  tip  indicated  that  the  rotor  would 
undergo  tip  stall  as  was  tested  and  described  in  "Measurement  of  the  Unsteady  Case  Wall 
Pressures."  (Gannon,  Hobson,  Shreeve,  2005) 
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3.  Comparison  of  Exit  Plane  CFD  Data  to  Experimental  Data 

The  CFD  data  produced  by  the  SWIFT  code  was  compared  to  experimental  data 
taken  from  downstream  of  the  rotor  for  the  equivalent  pressure  ratios.  The  following 
graphs  demonstrate  the  high  level  of  accuracy  of  the  CFD  for  the  90%  speed.  At  peak 
efficiency,  the  CFD  total  pressure  ratio  (Figure  43)  exhibited  nearly  the  same  form  and 
produced  data  very  close  to  that  of  the  experimental.  Both  Mach  number  calculations 
(Figure  44)  and  flow  angles  (Figure  45)  were  also  very  similar  and  all  three  graphs  are 
shown  below. 

However,  an  inaccuracy  was  introduced  during  the  production  of  the  CFD  data. 
In  the  experimental  model,  the  hub  did  not  rotate  downstream  of  the  rotor.  In  the  CFD 
grid,  a  rotating  hub  was  used  in  the  model  downstream  of  the  rotor,  producing 
discrepancies  in  the  flow  close  to  the  hub.  This  is  more  evident  in  the  flow  angle  graph 
shown  in  Figure  45.  Additional  graphs  for  all  settings  at  90%  speed  can  be  found  in 
appendix  M. 


Figure  43.  CFD  vs  experimental  data,  90%  speed,  peak  efficiency:  total  pressure  ratio. 
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Figure  44.  CFD  vs  experimental  data,  90%  speed,  peak  efficiency:  Mach  number. 
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Figure  45.  CFD  vs  experimental  data,  90%  speed,  peak  efficiency:  flow  angle. 
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4.  Comparison  of  CFD  Diffusion  Factor  Data  to  Experimental  Results 

The  diffusion  factor  results  were  graphed  and  compared  to  the  SWIFT  code 
results.  The  CFD  data  also  predicted  a  tip  stall.  Along  the  outer  half  of  the  blade  radius, 
the  CFD  and  experimental  data  exhibited  the  same  characteristics  and  obtained  very  close 
results  for  both  near  stall  and  peak  efficiency.  The  choke  data,  however,  showed  the 
largest  discrepancy.  In  addition,  in  the  inner  half  of  the  blade  radius,  the  experimental 
and  CFD  data  exhibited  different  characteristics,  although  the  magnitude  of  the  data  was 
similar.  This  also  can  be  attributed  to  the  rotating  hub  used  in  the  geometric  grid  in  the 
CFD  calculations.  Traveling  radially  inward,  the  experimental  data  decreased  towards 
zero,  while  the  CFD  data  decreased  by  the  same  magnitude  and  then  increased.  The 
results  are  shown  in  Figures  46-49. 


Hub-to-Tip  Radius  Ratio 


Figure  46.  70%  speed  diffusion  factor  comparison. 
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Figure  47.  80%  speed  diffusion  factor  comparison. 
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Figure  48.  90%  speed  diffusion  factor  comparison. 
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Figure  49.  100%  speed  diffusion  factor  comparison. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  present  study  set  up  the  instrumentation  and  established  a  baseline  for  future 
investigation  into  stall  and  surge  in  a  transonic  axial  compressor.  The  "pop  stall" 
phenomenon  experienced  by  jet  aircraft  ingesting  hot  gas  was  the  primary  purpose  of  the 
present  work. 

The  initial  probe  calibration  was  conducted  by  manually  rotating  the  probe  with 
an  electric  switch  controller.  Rotating  the  probe  to  a  specified  angle  using  the  manual 
controller  was  laborious  and  time  consuming. 

The  final  system  of  instrumentation  created  to  evaluate  flow  upstream  and 
downstream  of  the  compressor  reduced  the  amount  of  work  and  time  required.  The 
Matlab  control  system  using  USB  devices  worked  rapidly,  decreasing  the  time  required 
for  probe  movements.  The  pressure  transducers  provided  stable,  accurate  data  to  a  system 
capable  of  processing  the  data  as  it  was  taken  and  were  capable  of  being  calibrated  to  a 
lower  Mach  number.  The  USB  A/D  converter  allowed  rapid  communication  with  Matlab 
and  allowed  the  fast  processing  of  data  during  the  experiment.  In  addition,  the  USB  A/D 
converter  provided  feedback  for  the  probe  actuator  as  well  as  the  processing  of  pressure 
data  to  allow  feedback  for  to  null  yaw  the  probe.  However,  the  Matlab  programming 
could  be  optimized  as  the  system  slowed  occasionally  and  failed  when  the  probe  was 
placed  in  a  highly  swirling  flow. 

The  nature  of  the  diffusion  factor  calculated  was  found  to  be  close  to  what  was 
predicted  by  CFD  code  and  the  magnitude  was  within  the  expected  range.  In  addition, 
the  shape  of  the  diffusion  factor  data  shows  the  higher  loading  at  the  tip  and  validates  the 
prediction  of  a  tip  stall. 
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B.  RECOMMENDATIONS 

All  data  was  collected  with  a  single,  three-hole  probe.  To  calculate  the  diffusion 
factor,  data  from  both  upstream  and  downstream  was  necessary.  The  data  was  collected 
on  separate  days,  with  different  atmospheric  conditions.  Compressor  speed  and  mass 
flow  rates  were  repeated,  giving  close  approximate  reproductions  of  flow  conditions.  To 
ensure  less  error,  the  use  of  two  probes,  giving  simultaneous  data  upstream  and 
downstream  is  recommended. 

Data  downstream  of  the  rotor  was  collected  with  a  three-hole  probe  that  was  not 
sensitive  to  pitch.  CFD  models  have  shown  that  the  flow  at  the  exit  of  the  rotor  contains 
significant  radial  components,  as  well  as  axial  and  tangential  velocity  components.  A 
five-hole  probe  would  be  better  able  to  detennine  both  pitch  and  yaw  angles  and  give  a 
better  representation  of  the  downstream  flow. 

The  temperature  data  used  in  this  experiment  was  provided  by  thermocouples 
connected  to  the  HPVEE  system  run  from  a  separate  PC.  More  accurate  data  collection 
would  be  possible  with  a  thennocouple  connected  to  a  USB  device  controlled  with 
Matlab.  Temperature  measurements  would  then  be  collected  simultaneously  with  the 
pressure  data. 

The  diffusion  factor  calculations  require  the  relative  entrance  velocity  and  relative 
exit  velocity  of  the  flow  referenced  to  the  blade.  The  diffusion  factor  was  calculated 
from  holes  more  than  four  inches  upstream  and  two  inches  downstream  of  the  rotor. 
Gathering  data  closer  to  the  rotor  would  increase  the  accuracy  of  the  diffusion  factor 
calculations. 

The  Matlab  functions  used  to  control  the  probe  and  sample  data  could  be 
optimized.  A  better  logic  loop  for  rotating  to  an  angle  or  traversing  to  a  position  could 
reduce  the  time  required  for  each  experiment.  In  addition,  the  calling  of  calibration 
constants  required  a  large  amount  of  computer  memory  to  load  the  constants  in  the  MAT 
file  and  interpolate  the  angle  and  Mach  number  and  may  require  even  more  memory  for  a 
five-hole  probe. 
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APPENDIX  A:  CASE  WALL  ENGINEERING  DRAWING 
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Probe  Hole  Distances  from  Trailing  and  Leading  Edges  of  Rotor 


Hole  Number 

Hole  Location 

1 

1 1 1.5mm  (4.39")  Upstream  of  Rotor  Leading  Edge 

2 

9.91mm  (.39")  Upstream  of  Rotor  Leading  Edge 

3 

8.64mm  (.34")  Downstream  of  Rotor  Trailing  Edge 

4 

59.4mm  (2.34")  Downstream  of  Rotor  Trailing  Edge 
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Isentropic  effidency  T  otal  to  -total  pressure  ratio 


APPENDIX  B:  COMPLETE  STAGE  COMPRESSOR  MAP 
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APPENDIX  C:  COMPRESSOR  INTAKE  SYSTEM 
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APPENDIX  D:  SCREEN  SHOTS  OF  THE  HPYEE  PROGRAM 
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APPENDIX  E:  CALIBRATION  EQUATIONS  AND  SAMPLE 

CALCULATION 


In  the  calibration,  []  is  the  yaw  angle  taken  from  the  angle  of  the  probe  actuator, 
which  was  set  to  -20°  with  a  .  1°  tolerance. 


p  :=  -19.928 


The  pressures  P1-P3  were  taken  from  the  three-hole  probe. 


P1  :=  1 17780.4Pa 

P2  :=  1 17318. 3Pa 
P3  :=  101240. 6Pa 

P0j  and  Psj  were  taken  from  the  Stagnation  and  Static  Pressure  Probes. 

P0J-  :=  120177Pa 
PSJ-  :=  100783Pa 
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Tp  :=  297.5K 


Cp  :=  1005 


J 

k^K 


y  :=  1.41 


The  pseudo  dynamic  pressure,  q_probe  is  given  as: 

(p2  +  p3) 
q_probe  :=  P1 - - - 

q_probe  =  8.501  x  103Pa 


VJet  := 


m 

VJet  =  172.72  — 
s 


The  angle  flow  discriminators  give  the  direction  of  the  flow. 
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(Pl-P2) 

q_probe 


Cp  i  =  0.054 


(Pl-P3) 

q_probe 


Cp_2=  1.946 


(P2-P3) 

q_probe 


Cp_3  =  1-891 


The  total  and  static  pressure  coefficients  provide  the  measurable  data  for 
discerning  unknown  flows. 
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\  Probe  ’ 


(pl-p0j) 

q_probe 


At  Probe  -  -°-282 


Probe  • 


q_probe 
(Poj  “  Psj) 


As_Probe  0-438 


(y— i) 


MachNo  := 


Y 


fjJ 

vY-  lj 


MachNo  =  0.506 


76 


APPENDIX  F:  MATLAB  CALIBRATION  FILES 


function^  =  calibration_run(angle1 ,  angle2,  points,  Tp,  Tj) 
positions=linspace(angle1 ,  angle2,  points); 

for  i=1  :points 

lnput_Angle=positions(i); 

handles  =  guihandles(gcbo);  %  generate  handles  struct 

set(handles.lnput_Angle,  'String',  lnput_Angle); 

des_angle=lnput_Angle; 

rotate(des_angle); 

pause(2); 

OUTPUT_TO_EXCEL(Tp,  Tj); 
end 


function[cur_angle]  =  rotate(des_angle) 

%this  function  calls  the  data  from  the  GUI,  PMD_1608FS,  and  tells  the 

%ERB_24  to  rotate  the  probe  actuator.  The  desired  angle  (des_angle)  is  input  from  the 

GUI. 


%For  Loop  to  get  to  the  proper  angle 
for  i=1 :15; 

%Get  the  current  analog  angle  and  call  it  'h'. 
h=PMD_1 608FS_2; 

%Convert  the  analog  angle  into  a  real  angle. 
cur_angle=angle_calibration(h); 


handles  =  guihandles(gcbo);  %  generate  handles  struct 
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%Set  the  value  of  the  Cur_Angle  textbox  to  the  current  angle. 
set(handles.Cur_Angle,  'String',  cur_angle); 

%Assign  the  cur_angle  to  be  used  in  the  base  workspace. 
assignin('base','cur_angle',cur_angle) 

%Find  the  difference  between  the  current  angle  and  the  required  angle. 
diff=cur_angle-des_angle; 
if  abs(diff)<=.  1 
break; 
end 


if  cur_angle>=des_angle&&abs(diff)>=20 

ERB24_LEFT(6); 

pause(.l) 

elseif  cur_angle>=des_angle&&abs(diff)<=1 9.99999&&abs(diff)>=1 0 

ERB24_LEFT(3); 

pause(.l) 

elseif  cur_angle>=des_angle&&abs(diff)<=9.99999&&abs(diff)>=5 
ERB24_LEFT(1 .5); 


elseif  cur_angle>=des_angle&&abs(diff)<=4.99999&&abs(diff)>=3 
ERB24_LEFT(1) 

elseif  cur_angle>=des_angle&&abs(diff)<=2.99999&&abs(diff)>=1 
ERB24_LEFT(.1) 

elseif  cur_angle>=des_angle&&abs(diff)<=.999999&&abs(diff)>=.2 
ERB24_LEFT(.02) 

elseif  cur_angle>=des_angle&&abs(diff)<=.1 99999&&abs(diff)>=.1 5 
ERB24_LEFT(.01 ) 

elseif  cur_angle>=des_angle&&abs(diff)<=.1499999&&abs(diff)>=.1 
ERB24_LEFT(.003) 


end 


if  cur_angle<=des_angle&&abs(diff)>=20 

ERB24_RIGHT(6); 

pause(.l) 

elseif  cur_angle<=des_angle&&abs(diff)<=1 9.99999&&abs(diff)>=1 0 
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ERB24_RIGHT(3); 

pause(.l) 

elseif  cur_angle<=des_angle&&abs(diff)<=9.99999&&abs(diff)>=5 
ERB24_RIGHT(1 .5); 


elseif  cur_angle<=des_angle&&abs(diff)<=4.99999&&abs(diff)>=3 
ERB24_RIGHT(1) 

elseif  cur_angle<=des_angle&&abs(diff)<=2.99999&&abs(diff)>=1 
ERB24_RIGHT(.1) 

elseif  cur_angle<=des_angle&&abs(diff)<=.999999&&abs(diff)>=.2 
ERB24_RIGHT(.02) 

elseif  cur_angle<=des_angle&&abs(diff)<=.1 99999&&abs(diff)>=.1 5 
ERB24_RIGHT(.01) 

elseif  cur_angle<=des_angle&&abs(diff)<=.149999&&abs(diff)>=.1 
ERB24_RIGHT(.003) 


function[current_posit]  =  angle_calibration(analog_posit) 

%insert  data  into  x  and  y  to  get  the  calibration  curve  for  rotation. 
%the  x  data  used  for  the  angle  was  taken  from  digital  inclinometer 
%the  y  data  was  read  from  the  digital  pmdl  608fs_2 


%  First  calibration  prior  to  connecting  all  info  through  Amphenol  cord. 
%(4/1 1/05) 


%x=[-50.3  -40.3  -29.9  -19.8  -10.2  0  10  20.3  30.1  40.2  50.1]; 

%y=[1 .093  1 .350  1.641  1.906  2.160  2.432  2.712  3  3.263  3.530  3.79]; 

%Second  Calibration  Data  after  connecting  through  Amphenol  cord. 

%(4/1 3/05) 

x=[-49.8  -40  -29.9  -19.95  -9.8  .2  10.4  20  30.2  40.3  50.4]; 
x=x+.4; 

y=[1. 0978  1.3665  1.6335  1.9022  2.1701  2.4381  2.7086  2.9846  3.2541  3.5261  3.8]; 

%create  a  square  polynomial  to  fit  the  data  to. 
p=polyfit(y,x,2); 

%evaluate  the  polynomial  at  the  analog  position  and  output  the  position  in 
%degrees 

current_posit=polyval(p,analog_posit); 
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function[analog_pressure]=PMD_1608FS_3 

%This  function  is  called  to  take  an  average  of  channel  one's  output. 

%lt  will  be  used  for  feedback  on  the  pressure  of  the  three  hole  probe. 

clear  all 

global  analog_pressure 


openDAQ=daqfind; 
for  i=1  :length(openDAQ), 
stop(openDAQ(i)); 
end 


%stop(ai); 

%delete(ai); 

%clear  all 

ai  =  analoginput('mcc',0); 

%create  channel  3-7  (called  2-6) 
addchannel(ai,2:6); 

set(ai,'SampleRate',500); 

ActualAIRate  =  ai.SampleRate; 
set(ai,'SamplesPerTrigger',30) 

%the  t  is  left  in  for  plotting  and  troubleshooting 
start(ai) 

[a,t]  =  getdata(ai); 
start(ai) 

[b,t]  =  getdata(ai); 
start(ai) 

[c,t]  =  getdata(ai); 
start(ai) 

[d,t]  =  getdata(ai); 

%plot(a,t) 

%take  an  average  of  channel  each  output. 

analog_pressure=(sum(a)/30+sum(b)/30+sum(C)/30+sum(d)/30)/4; 


%stop  and  delete  the  ai 

stop(ai); 

delete(ai); 
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function[Pressures]  =  pressure_calibration 

%insert  data  into  x  and  y  to  get  the  calibration  curve  for  rotation. 
%the  x  data  used  for  the  angle  was  taken  from  digital  inclinometer 
%the  y  data  was  read  from  the  digital  pmdl  608fs_2 


%  First  calibration. 

%(5/3/05) 

%pressure  in  inches  of  mercury  converted  to  Pascals 

global  Pressures 
global  analog_pressure 

analog  _pressure=PMD_1 608FS_3; 

x=[0  2  4  6  8  10  12  14  1 6]'*3386.38+[29.795]*3386.38; 


%Values  for  upper  pressure  connection. 


%y(:,1)=[3. 50302383 
%3. 339436942 
%3. 172269271 
%3.01 0462603 
%2. 844668243 
%2. 679341 828 
%2.51 9305206 
%2. 354955367 
%2.1 91 195544 
%]; 

%y(:,2)=[3. 50691 9966 
%3.341 28837 
%3. 171496147 
%3. 007471 834 
%2. 83968363 
%2.671 834389 
%2. 509325806 
%2. 342646423 
%2.1 7691 31 
%]; 

%y(:,3)=[3.5041 93688 
%3. 342549783 
%3.1 760331 63 
%3.01 58541 24 
%2. 850883752 
%2. 686930648 
%2. 52751 456 
%2. 364995804 
%2. 202263422 
%]; 

%y(:,4)=[3.51 2240279 
%3. 346832481 
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%3. 17700974 
%3.01 3300781 
%2. 845888965 
%2.6781 71 969 
%2.51 57651 13 
%2.3491 06076 
%2. 183474479 
%]; 


%Values  for  Upper  pressure  connection  5_5_05 


y(:,1)=[3.5036 

3.3389 

3.1709 

3.0046 

2.8464 

2.6802 

2.5208 

2.3544 

2.1891 

]; 

y(:,2)=[3.5082 

3.341 

3.1703 

3.0016 

2.8412 

2.6727 

2.5107 

2.3419 

2.1742 

]; 

y(:,3)=[3.5048 

3.3418 

3.1744 

3.0098 

2.8527 

2.688 

2.5292 

2.3641 

2.2001 

]; 

y(:,4)=[3.5 

3.3348 

3.1685 

3.0069 

2.8396 

2.6759 

2.5129 

2.3494 

2.1863 

]; 


y(:,5)=[3.5129 

3.3457 

3.1751 
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3.0065 

2.8469 

2.6785 

2.5167 

2.3473 

2.1803 

]; 


for  ii=1 :5 

%create  a  polynomial  to  fit  the  data  to.  Data  was  nearly  linear. 
p(ii,:)=polyfit(y(:,ii),x,1); 

%evaluate  the  polynomial  at  the  analog  position  and  output  the  position  in 
%degrees 

Pressures(ii)=polyval(p(ii,:),analog_pressure(ii)); 

end 


function[]=OUTPUT_TO_EXCEL(Tp,  Tj) 

%this  function  gets  the  data  from  the  devices  using  PMD_1608FS_2  and 
%PMD_1608FS_3  and  sends  the  data  to  ExcelWrite  to  be  written  to  an  Excel 
%Spreadsheet. 


global  Pressures 


%Freddy.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  headers  will  be  created, 
load  freddy.txt 

%restart.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  data  will  not  be  overwritten  if  the  run  is  restarted, 
load  restart.txt 

%if  it  is  the  first  run,  or  the  run  has  been  restarted,  this  loop  will 
%place  the  header  in  and  then  write  the  first  line  of  data, 
if  freddy==1 ; 


%if  first  line,  put  the  header  in.  if  restart,  add  the  number  of 
%lines  so  that  data  will  not  be  overwritten. 
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a=freddy+restart; 


%header 

m={'angle',  'PI',  'P2',  'P3',  'Poj',  'Psj',  'Tp',  'Tj',  'q_probe',  'VJet',  'Cp1_Probe',... 
’Cp2  Probe’,  'Cp3_Probe',  'At_Probe',  'As_Probe',  'MachJSIo'}; 


%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

fred  =  [’ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls 

. ['  num2str(a) num2str(2) '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 
eval(fred); 


%change  the  row  number  it  will  be  the  next  available. 


%Call  the  calculation  function 
calculate(Tp,  Tj); 

%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

%fred  =  ['ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls 

0,'... 

%  . ['  num2str(a) ','  num2str(2) '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 

%eval(fred); 

%Change  the  value  of  freddy  to  the  next  row  number. 

%freddy=freddy+1 ; 

%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared. 

%save  freddy.txt  freddy  -ascii 


%lf  the  run  is  not  the  first,  or  a  reset,  just  output  dats. 
elseif  freddy  ~=1 

a=freddy+restart; 


%Call  the  calculation  function 
calculate(Tp,  Tj); 
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%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

%fred  =  ['ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls' 

%  . ['  num2str(a) num2str(2)  '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 

%eval(fred); 

%Change  the  value  of  freddy  to  next  row  number. 

%freddy=freddy+1 ; 

%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared. 

%save  freddy.txt  freddy  -ascii 
end 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%%% 

%SUBFUNCTIONS 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/ 

/o  /o  /o 


function[m]  =  calculate(Tp,  Tj) 

%this  subfunction  will  get  all  required  data  three  times  and  average 
%the  three  sets  of  values  before  outputting  the  final  answer. 


%%%%%%%%%%%%CONSTANTS%%%%%%%%%%%%%%%%%% 


%Cp  of  air 
Cp_air=1005; 


%gamma  of  air 

gam_air=1 .41 ;  %gamma  of  air 


%gravitational  constant  (m/secA2) 
g=9.81 ;  %gravitational  constant  (m/secA2) 


%density  of  Mercury  (kg/mA3)  (calibration  pressure  in  inches  of  mercury) 
rho_Hg= 13550; 


%loop  for  data 
for  i=1 :3 

%Get  the  data  from  the  PMD  1608FS 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o 

%Get  the  current  angle. 
h=PMD_1608FS_2; 
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%Convert  the  angle  to  degrees. 
cur_angle(i)=angle_calibration(h); 

%get  the  pressure  variable  info. 

Pressures=pressure_calibration; 
%Pressures=PMD_1 608FS_3; 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%Cobra  Probe  pressures 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o 

%Left  Hole  Pressure 
P2=Pressures(1); 


%Center  Hole  Pressure 
P1=Pressures(2); 


%Right  Hole  Pressure 
P3=Pressures(3); 

%Total  Pressure  of  jet 
Poj=Pressures(5); 

%Static  Pressure  of  jet 
Psj=Pressures(4); 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%Temperature  in  pipe 

%Tp=is  entered  manually  in  the  GUI; 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%%%%%%%%%%%%CALCULATIONS%%%%%%%%%%%%%%%%%%%% 


q_probe=P1-(P2+P3)/2; 

VJet(i)=sqrt(1-(Psj/Poj)A((gam_air-1)/gam_air))*sqrt(2*Cp_air*Tp); 

Cp1_Probe(i)=(P1  -P2)/q_probe; 

Cp2_Probe(i)=(P1  -P3)/q_probe; 

Cp3_Probe(i)=(P2-P3)/q_probe; 

At_Probe(i)=(P1  -Poj)/q_probe; 
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As_Probe(i)=q_probe/(Poj-Psj); 

%Mach  Number 

Mach_No(i)=sqrt(((Poj/Psj)A((gam_air-1  )/gam_air)-1  )*(2/(gam_air-1 ))); 
end 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%Average  the  values 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


cur_angle=(cur_angle(1)+cur_angle(2)+cur_angle(3))/3; 

Cp1_Probe=(Cp1_Probe(1)+Cp1_Probe(2)+Cp1_Probe(3))/3; 

Cp2_Probe=(Cp2_Probe(1)+Cp2_Probe(2)+Cp2_Probe(3))/3; 

Cp3_Probe=(Cp3_Probe(1)+Cp3_Probe(2)+Cp3_Probe(3))/3; 

VJet=(VJet(1)+VJet(2)+VJet(3))/3; 

At_Probe=  (At_Probe(1  )+At_Probe(2)+At_Probe(3))/3; 

As_Probe=(As_Probe(1)+As_Probe(2)+As_Probe(3))/3; 

Mach_No=(Mach_No(1)+Mach_No(2)+Mach_No(3))/3; 


Tj=Tj; 


%Create  the  data  array  to  be  written. 

m=[cur_angle,  Pressures(3),Pressures(2),Pressures(1),  Pressures(5),  Pressures(4),  Tp, 
Tj,  q_probe,  VJet,  Cp1_Probe,... 

Cp2_Probe,  Cp3_Probe,  At_Probe,  As_Probe,  Mach_No]; 

%Freddy.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  headers  will  be  created, 
load  freddy.txt 

%restart.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  data  will  not  be  overwritten  if  the  run  is  restarted, 
load  restart.txt 

a=freddy+1+restart; 

%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

fred  =  ['ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls . 

. ['  num2str(a)  7  num2str(2) '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 
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eval(fred); 


%Change  the  value  of  freddy  to  next  row  number. 
freddy=freddy+1 ; 

%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared, 
save  freddy.txt  freddy  -ascii 
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APPENDIX  G:  SCANIVALVE  CALIBRATION  GRAPHS 
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APPENDIX  H:  PX-138  CALIBRATION  GRAPHS 
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APPENDIX  I:  MATLAB  INTERPOLATION  FILES 


%  m-file  to  call  the  3  hole  probe  output  file 
clear  all 


%  Type  of  file  to  get  data  from 
%file_type  =  ' e';  %  Call  from  excel 

file_type  =  'm';  %  Call  from  text  file 

%  Measured  values,  can  be  changed  to  a  function  later  on 

pl_m  =  [93692.44746 

93648 . 41833 

93540 . 03893 

93419 . 38217 

93308 . 88599 

93144 . 62346 

93158 . 17088 

93047 . 25134 

93104 . 40454 

93128 . 95925 

93103.98118 

93125.99575 

92763.60213 

92046 . 01196 

91071 .86744 

] ; 

p2_m  =  [95925.86792 
95915.51443 
95922 . 96894 
95927 . 52448 
95932 . 08001 
95862 . 50458 
95918 . 82755 
95897 . 70643 
95905.98922 
95862 . 91872 
95831 . 44412 
95680 . 69734 
94890 . 93334 
93210 . 3553 
91121 . 43572 
] ; 

p3_m  =  [93731.43686 

93648 . 22584 

93532 .23472 

93429.69185 

93297 . 73094 

93106.09345 

93086.76159 

93009.85444 

93100 . 6301 

93164 . 089 

93188 .46395 
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93155.26359 
92857 . 30091 
92044 . 94283 
91100 . 62385 

] ; 

poj_m  =  [95852.60365 

95850 . 08642 

95877 . 35637 

95901 .27002 

95863.51163 

95820 . 71878 

95880 . 71268 

95850 . 9255 

95881 . 97129 

95873.58053 

95845 . 47151 

95879.03452 

95921 . 40784 

95846 . 73012 

95884 .06898 

] ; 


%pl_m 
%p2_m 
%p3_m 
%po j_m 


102794 . 6236; 
107286.9565; 
105166.2278; 
107363.8443; 


tic 

% [Pitch_m, Yaw_m, Mach_m, ps_m, po_m] 
three_hole_probe_output (f ile_type, pl_m, p2_m, p3_m, po j_m) 


[Pitch_m_2 , Yaw_m_2 , Mach_m_2 , ps_m_2 , po_m_2 ] 
three_hole_probe_output ( f ile_type, pl_m, p2_m, p3_m, p2_m) 


toe 


%  M-function  file  to  obtain  the  flow  direction  and  mach  number 
based  on  the  input  data  from  a  3-hole  probe 
"6 

%  To  assume  ZERO  PITCH,  just  pass  p2_m  to  poj_m  so  that  they  are 
the  same.  This  will  stop  the  program 

%  searching  along  this  dimension. 

'O 

%  To  load  data  from  the  excel  spreadsheets  they  must  be  moved 
into  the  same  directory,  the  constants  file 
%  will  then  also  be  updated 
%  file_type  =  ' e'  (excek)  or  'm'  (matlab) 


function  [Pitch_m, Yaw_m, Mach_m, ps_m, po_m] 

three_hole_probe_output (f ile_type, pl_m, p2_m, p3_m, po j_m) 
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%  If  only  one  input  point  is  needed  two  points  are  required  to 
make  griddata3  work 

if  length (pl_m) <2 
point_no  =  1; 

pl_m  =  [pl_m  pl_m]  ;  p2_m  =  [p2_m  p2_m]  ;  p3_m  =  [p3_m 

p3_m] ;  poj_m  =  [poj_m  poj_m] ; 
else 

point_no  =  0; 

end 


%  Input  matricies  or  vectors  are  reshaped  to 
shape_pl_m  =  size (pl_m) ; 

reshape (pl_m, 1, shape_pl_m (1) *shape_pl_m (2) ) ; 

shape_p2_m  =  size (p2_m) ; 

reshape (p2_m, 1, shape_p2_m (1) *shape_p2_m (2) ) ; 

shape_p3_m  =  size (p3_m) ; 

reshape (p3_m, 1 , shape_p3_m ( 1 ) *shape_p3_m (2 ) ) ; 

shape_poj_m  =  size(poj_m); 

reshape (po j_m, 1 , shape_po j_m ( 1 ) *shape_po j_m (2 ) ) ; 


column  vectors 
pl_m 

p2_m 

p3_m 

po  j_m 


%  Physical  constants 
Gam_gas  =  1.41; 


%Measured  constants 

q_yaw_m  =  pl_m- ( (p2_m+p3_m) / 2)  ; 

q_pitch_m  =  poj_m- ( (p2_m+p3_m+pl_m) /3) ; 


Cpl_m 

Cp2_m 

Cp3_m 

Cp4_m 

Cp5_m 

Cp6_m 


(pl_m-p2_m) . /q_yaw_m; 

(pl_m-p3_m) . /q_yaw_m; 

(p3_m-p2_m) ,/q_yaw_m;  %  Side  holes 

(poj_m-p2_m) . /q_pitch_m; 

(poj_m-p3_m) . /q_pitch_m; 

(poj_m-pl_m) ,/q_pitch_m;  %  Centre  hole 


if  file_type  ==  ' e' 

%  Coefficients  of  pressure  are  read  in  but  are  calculated 
later  from  the  raw  data. 

%Cpl  =  Load_constant ( 12 ) 

%Cp2  =  Load_constant (13) 

%Cp3  =  Load_constant (14) 

Yaw  =  Load_constant (2 ) ; 

%  Pressure  are  read  in 
pi  =  Load_constant  ( 4 ) ; 
p2  =  Load_constant ( 3 ) ; 
p3  =  Load_constant ( 5 ) ; 
poj  =  Load_constant  (  6 ) ; 

%  At,  As  and  Mach  number  are  loaded  in 
At  =  Load_constant (15) ; 

As  =  Load_constant ( 1 6 )  ; 

Mach  =  Load_constant ( 17 )  ; 
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%  Pitches  are  constructed 
Pitch  =  Load_constant ( 18  )  ; 


save  3_hole_constants  Yaw  pi  p2  p3  poj  At  As  Mach  Pitch 

else 

load  3_hole_constants  Yaw  pi  p2  p3  poj  At  As  Mach  Pitch; 
end  %  if  file_type  ==  'excel' 

%  q  coefficients 

q_yaw  =  pi- (  (p2+p3) /2 )  ;  %  Same  one  as  is  used  in  the 

calibration 

q_pitch  =  po j- ( (p2+p3+pl) /3) ; 

%At  =  (pi  -  poj)  ./q_yaw  %  There  is  some  difference  between  the 
spreadsheet  and  these  values 

%  Pitch  presssure  coefficients  are  now  defined  for  a  '4'  hole 

probe 

Cpl  =  (pl-p2) . /q_yaw; 

Cp2  =  (pl-p3) . /q_yaw; 

Cp3  =  (p3-p2) . /q_yaw;  %  Side  holes 

Cp4  =  (po j-p2) . /q_pitch; 

Cp5  =  (po j-p3) . /q_pitch; 

Cp6  =  (poj-pl) ./q_pitch;  %  Centre  hole 

"6 


%  Initial  estimate  of  Mach  number 

ps_m  =  min ( [pl_m;  p2_m;  p3_m] ) ;  %ps_m  =  min ( [pl_m;  p2_m; 

p3_m;  poj_m] ) ; 

po_m  =  max ( [pl_m;  p2_m;  p3_m] ) ;  %po_m  =  max ( [pl_m;  p2_m; 

p3_m;  poj_m] ) ; 

Mach_m  =  sqrt ( (2/ (Gam_gas-1) ) * ( (ps_m. /po_m) . A ( (1- 

Gam_gas ) /Gam_gas ) -1 ) ) ; 

Pitch_m  =  0*ones (size (pl_m) ) ;  %  Initial  estimate  of  pitch  in 

degrees 

Yaw_m  =  0*ones (size (pl_m) ) ;  %  Initial  estimates  of  pitch 


for 

Mach_m ( i ) ) 
Cp3_m ( i ) ) 
Cp  6_m ( i ) ) 


i  =  1 : length (Mach_m) 
if  (max (max (Mach) )  < 

I  •  •  • 

(max (max (Cp3 ) ) 

I  -  -  ■ 

(max (max (Cp6 ) ) 


Mach_m (i) ) 

<  Cp3_m ( i ) ) 

<  Cp  6_m ( i ) ) 


(min (min (Mach) )  > 

(min (min (Cp3 ) )  > 

(min (min (Cp6 ) )  > 


disp (' Warning,  data  may  be  outside  calibration  range') 
Yaw_m(i)  =  0/0; 

end  %  if  (max (max (Mach) )  <  max (Mach_m) ) 

end 
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temp  =  f ind ( ( isnan ( Yaw_m) ) ~=1 )  ;  %  Data  that  is  within  calibration 

range 


%  Initial  estimate  of  yaw 

Yaw_m(temp)  = 

griddata3 (Mach, Pitch, Cp3, Yaw, Mach_m (temp) , Pitch_m (temp) , Cp3_m (temp) )  ; 

%  Variables  on  the  tolerance 

tol  =  0.25;  Pitch_m_old  =  Pitch_m+2 *tol ;  Yaw_m_old  = 

Yaw_m+2*tol ; 

count  =  0; 

while  (max ( [max (abs ( Y a w_m- Y a w_m_o 1 d ) )  max (abs (Pitch_m- 

Pitch_m_old) ) ] )  >  tol)  &  count  <  10 

%  Stored  variables  for  the  tolerance 

Pitch_m_old  =  Pitch_m;  Yaw_m_old  =  Yaw_m; 

count  =  count+1; 

if  length  (intersect  (poj_m,  p2_m)  )  ==  0  %  is  this  is  true  then 

it  is  assumed  that  the  pitch  is  zero 

%  Pitch  is  calculated  using  the  centre  hole  but  in  case 
this  is  outside  the  range  the  other  two  are  used 

Pitch_m (temp)  = 

griddata3 (Yaw, Cp6, Mach, Pitch, Yaw_m (temp) , Cp6_m (temp) , Mach_m (temp) ) 

if  isnan (Pitch_m (temp) ) 

Pitch_m (temp)  = 

griddata3 (Yaw, Cp5, Mach, Pitch, Yaw_m (temp) , Cp5_m (temp) , Mach_m (temp) ) ; 
end 

if  isnan (Pitch_m (temp) ) 

Pitch_m (temp)  = 

griddata3 (Yaw, Cp4 , Mach, Pitch, Yaw_m (temp) , Cp4_m (temp) , Mach_m (temp) ) ; 
end 

end 

%  Yaw  is  corrected  to  the  right  pitch 

Yaw_m(temp)  = 

griddata3 (Mach, Pitch, Cp3, Yaw, Mach_m (temp) , Pitch_m (temp) , Cp3_m (temp) ) ; 

%  The  static  and  stagnation  pressure  constants  are  found 

At_m(temp)  = 

griddata3 (Mach, Pitch, Yaw, At , Mach_m (temp) , Pitch_m (temp) , Yaw_m (temp) )  ; 

As_m(temp)  = 

griddata3 (Mach, Pitch, Yaw, As, Mach_m (temp) , Pitch_m (temp) , Yaw_m (temp) ) ; 

%  The  total  pressure  from  the  constants  is  calculated  even 
though  for  the  3-hole  probe  it  is  directly  measured 

po_m(temp)  =  p2_m(temp)  -  At_m (temp) . *q_yaw_m (temp) ; 

%  The  static  pressure  is  calculated 

ps_m(temp)  =  po_m(temp)  -  q_yaw_m (temp) . /As_m (temp) ; 

%  The  mach  number  is  updated 
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Mach_m (temp)  =  sqrt ( (2/ (Gam_gas- 

1 ) ) * ( (ps_m (temp) . /po_m (temp) ) . A ( ( 1-Gam_gas ) /Gam_gas ) -1 ) ) ; 

end  %  while  abs (Pitch_m-Pitch_m_old) >tol 

%  Points  are  reshaped  to  the  original  input  matrix 
if  point_no  ==  1 


Pitch_m 

=  Pitch_m ( 1 ) ; 

Yaw_m 

=  Yaw_m ( 1 ) ; 

Mach_m 

=  Mach_m ( 1 ) ; 

ps_m 

=  ps_m (1) ; 

po_m 

=  po_m (1) ; 

Pitch_m 

=  reshape (Pitch_m, shape_pl_m) ; 

Yaw_m 

=  reshape (Yaw_m, shape_pl_m) ; 

Mach_m 

=  reshape (Mach_m, shape_pl_m) ; 

ps_m 

=  reshape (ps_m, shape_pl_m) ; 

po_m 

=  reshape (po_m, shape_pl_m) ; 

end  %  if  point_no  ==  1 


if  0 

%  Cpl  and  Cp2  are  plotted 
figure  (1);  close;  figure  (1); 
plot ( Yaw, Cpl , Yaw, Cp2 ) 

temp  =  axis;  axis ( [-30  30  temp (3)  temp (4)]); 

%  Cp3  is  plotted 

figure  (2);  close;  figure  (2); 

plot (Yaw, Cp3 ) 

temp  =  axis;  axis ( [-30  30  temp (3)  temp (4)]); 

%  Cp4  and  Cp6  is  plotted 

figure  (3);  close;  figure  (3); 

plot (Pitch ' , Cp4 ' , Pitch ' , Cp5 ' ) 

temp  =  axis;  axis ( [-30  30  temp (3)  temp (4)]); 


%  Cp6  is  plotted 

figure  (4);  close;  figure  (4); 

plot (Pitch ' , Cp6 ' ) 

temp  =  axis;  axis ( [-30  30  temp (3)  temp (4)]); 


%  Surface  plot  of  the  data  is  plotted 
figure  (5);  close;  figure  (5); 

TRI  =  delaunay (Yaw, Pitch) ; 
trisurf (TRI, Yaw, Pitch, Cp 3) 
hold  on 

%plot3 (Yaw_m, Pitch_m, Cp3_m,  ' or ' ) 
end  %  if  0 


%  m-function  file  to  pull  in  certain  data  sets 
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function  [A]  =  Load_constant (column) 


%  Mach  0.06  Nominal 
A  ( : ,  1 ) 

xlsreadcells ( ' 5_1 0_0  5_4  Odeg_Pitch . xls ' ,  ' Sheet  1 ' ,  [ column] ,  [110:122] )  ; 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls',  ' Sheet  1 ' ,  [ column] ,  [1:13]); 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_2  0deg_Pitch .xls',  ' Sheet  1 '  ,  [ column] ,  [1:13]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls',  ' Sheet  1 '  ,  [ column] ,  [1:13]); 

%  10  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' ,  ' Sheet 1 ' ,  [ column] ,  [111:123] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_10_05_10deg_Pitch_forward. xls', ' Sheetl ' , [column] , [104:1 
16] ) ;  %  -10  Degree  pitch 

%  Mach  0.15  Nominal 
A  (  : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_4  0deg_Pitch .xls',  'Sheetl',  [ column] ,  [92:104]); 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls', 'Sheetl', [ column] , [15:27]); 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_20deg_Pitch .xls', 'Sheetl', [ column] , [15:27]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls', 'Sheetl', [ column] , [16:28]); 

%  10  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' , ' Sheetl ' , [ column] , [94:106] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_10deg_Pitch_f orward .xls',  'Sheetl',  [ column]  ,  [87:99 
] ) ;  %  -10  Degree  pitch 

%  Mach  0.22  Nominal 
A ( : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_4  0deg_Pitch .xls',  'Sheetl',  [ column] ,  [76:88]); 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls', 'Sheetl', [ column] , [29:41]); 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_20deg_Pitch .xls', 'Sheetl', [ column] , [29:41]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls', 'Sheetl', [ column] , [30:42]); 

%  10  Degree  pitch 
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A  (  :  , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' , ' Sheetl ' , [column] , [77:89]); 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_10_05_10deg_Pitch_forward. xls', 'Sheetl', [column] , [69:81 
] ) ;  %  -10  Degree  pitch 

%  Mach  0 . 3  Nominal 
A  (  : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_4 Odeg_Pitch .xls',  'Sheetl',  [ column] ,  [3:15]); 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls',  'Sheetl',  [ column] ,  [43:55]); 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_20deg_Pitch .xls', 'Sheetl', [ column] , [43:55]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls', 'Sheetl', [ column] , [45:57]); 

%  10  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' , ' Sheetl ' , [ column] , [1:13] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_10_05_10deg_Pitch_forward. xls', 'Sheetl', [column] , [1:13] 

) ;  %  -10  Degree  pitch 

%  Mach  0 . 5  Nominal 
A ( : , end+1 ) 

xlsreadcells ( ' 5_10_05_4  0deg_Pitch .xls',  'Sheetl',  [ column] ,  [22:34] ) ; 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls', 'Sheetl', [ column] , [60:72]); 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_20deg_Pitch .xls', 'Sheetl', [ column] , [60:72]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls', 'Sheetl', [ column] , [64:76]); 

%  10  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' ,  ' Sheetl ' ,  [column] ,  [21:33] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_10_05_10deg_Pitch_forward. xls', 'Sheetl', [column] , [17:29 
] ) ;  %  -10  Degree  pitch 

%  Mach  0.7  Nominal 
A  (  : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_4  0deg_Pitch .xls',  'Sheetl',  [ column] ,  [39:51]); 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls', 'Sheetl', [ column] , [77:89]); 

%  30  Degree  pitch 
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A  (  :  , end+1 ) 

xlsreadcells ( ' 5_9_05_2  0deg_Pitch . xls ' ,  ' Sheet  1 ' ,  [ column] ,  [77:89]); 

%  20  Degree  pitch 
A  (  : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls',  ' Sheet  1 ' ,  [ column] ,  [82:94]); 

%  10  Degree  pitch 
A  (  :  , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' , ' Sheet 1 ' , [column] , [40:52] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_10deg_Pitch_f orward .xls',  ' Sheet  1 ' ,  [ column] ,  [34:46 
] ) ;  %  -10  Degree  pitch 

%  Mach  0.85  Nominal 
A ( : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_4  0deg_Pitch .xls',  ' Sheet  1 ' ,  [ column] ,  [59:71]); 

%  40  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_30deg_Pitch .xls', ' Sheetl ' , [column] , [93:105] ) ; 

%  30  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_20deg_Pitch .xls', 'Sheetl', [column] , [94:106]); 

%  20  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_9_05_10deg_Pitch .xls', 'Sheetl', [ column] , [101:113] ) ; 

%  10  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_5_05_0deg_Pitch . xls ' , ' Sheetl ' , [column] , [59:71] ) ; 

%  0  Degree  pitch 
A ( : , end+1 ) 

xlsreadcells ( ' 5_1 0_05_10deg_Pitch_f orward .xls', 'Sheetl', [ column] , [52:64 
] ) ;  %  -10  Degree  pitch 

%A  =  reshape (A, [ size (A, 1 ) *size (A, 2 )  1]); 
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APPENDIX  J:  UPSTREAM  PROBE  DATA:  STAGNATION 

PRESSURE 


70  %  Rotor  Speed 


80%  Rotor  Speed 


Probe  Position  4m} 
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90%  Rotor  Speed 


100%  Rotor  Speed 
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APPENDIX  K:  UPSTREAM  PROBE  DATA:  YAW  ANGLE 

DISTRIBUTION 


70  %  Rotor  Speed 


80  %  Rotor  Speed 
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90  %  Rotor  Speed 


100  %  Rotor  Speed 
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APPENDIX  L:  DOWNSTREAM  PROBE  DATA:  STAGNATION 

PRESSURE  RATIO 
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APPENDIX  M:  CFD  VS  EXPERIMENTAL  DATA:  90%  ROTOR 

SPEED 

Near  Stall  Rotor  Exit  Flow  Angle 


Near  Stall  Rotor  Exit  Mach  Number 
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Near  Stall  Rotor  Exit  Total  Pressure  Ratio 
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Peak  Efficiency  Rotor  Exit  Flow  Angle 


114 


Peak  Efficiency  Rotor  Exit  Mach  Number 


Peak  Efficiency  Rotor  Exit  Total  Pressure  Ratio 
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Choke  Rotor  Exit  Flow  Angle 


Choke  Rotor  Exit  Mach  Number 
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Choke  Rotor  Exit  Total  Pressure  Ratio 
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APPENDIX  N:  MATLAB  EXPERIMENTAL  DATA  COLLECTION 

FILES 


function[]=OUTPUT_TO_EXCEL(Tp,  Tj) 

%this  function  gets  the  data  from  the  devices  using  PMD_1608FS_2  and 
%PMD_1608FS_3  and  sends  the  data  to  ExcelWrite  to  be  written  to  an  Excel 
%Spreadsheet. 


global  Pressures 


%Freddy.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  headers  will  be  created, 
load  freddy.txt 


%restart.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  data  will  not  be  overwritten  if  the  run  is  restarted, 
load  restart.txt 


%if  it  is  the  first  run,  or  the  run  has  been  restarted,  this  loop  will 
%place  the  header  in  and  then  write  the  first  line  of  data, 
if  freddy==1 ; 


%if  first  line,  put  the  header  in.  if  restart,  add  the  number  of 
%lines  so  that  data  will  not  be  overwritten. 
a=freddy+restart; 


%header 

header={'probe_pos', 'angle',  'PI',  V2\  'P3',  'Poj',  'Psj',  'Tp',  'Tj',  'q_probe',  'VJet', 
'Cp1_Probe',... 

'Cp2_Probe',  'Cp3_Probe',  'At_Probe',  'As_Probe', 

'MachJ9o','Pitch_m','Yaw_m','Mach_m','ps_m','po_m', 
'Pitch_m_2','Yaw_m_2','Mach_m_2','ps_m_2','po_m_2','flow_angle'}; 


%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

fred  =  ['ExcelWrite('  ""  'C:\Documents  and 

Settings\user\Desktop\Transonic\testing1  .xls . ,  0,'... 

. ['  num2str(a)  ','  num2str(2) '] . ,  header)’]; 
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%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 
eval(fred); 


%change  the  row  number  it  will  be  the  next  available. 


%Call  the  calculation  function 
calculate(Tp,  Tj); 


%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

%fred  =  ['ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls . 

o;... 

%  . ['  num2str(a) num2str(2) '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 

%eval(fred); 

%Change  the  value  of  freddy  to  the  next  row  number. 

%freddy=freddy+1 ; 

%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared. 

%save  freddy.txt  freddy  -ascii 


%lf  the  run  is  not  the  first,  or  a  reset,  just  output  dats. 
elseif  freddy  ~=1 

a=freddy+restart; 

%Call  the  calculation  function 
calculate(Tp,  Tj); 

%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

%fred  =  [’ExcelWrite( . C:\Documents  and  Settings\user\Desktop\ERB24\testing1.xls’ 

" ',  0,'... 

%  . ['  num2str(a)  7  num2str(2)  '] . ,  m)']; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 
%send  it  to  an  EXCEL  Spreadsheet. 

%eval(fred); 

%Change  the  value  of  freddy  to  next  row  number. 

%freddy=freddy+1 ; 
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%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared. 
%save  freddy.txt  freddy  -ascii 
end 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%%% 

%SUBFUNCTIONS 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/ 

/o  /o  /o 


function[header]  =  calculate(Tp,  Tj) 

%this  subfunction  will  get  all  required  data  three  times  and  average 
%the  three  sets  of  values  before  outputting  the  final  answer. 


%%%%%%%%%%%%CONSTANTS%%%%%%%%%%%%%%%%%% 


%Cp  of  air 
Cp_air=1005; 


%gamma  of  air 

gam_air=1 .41 ;  %gamma  of  air 


%gravitational  constant  (m/secA2) 
g=9.81 ;  %gravitational  constant  (m/secA2) 


%density  of  Mercury  (kg/mA3)  (calibration  pressure  in  inches  of  mercury) 
rho_Hg= 13550; 

for  k=1 :2 


%loop  for  data 
for  i=1 :3 


%Get  the  data  from  the  PMD  1 608FS's 


O/  O /  O /  O/  O/  O/  O/  Of  o /  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o 

%Get  the  current  angle. 
h=PMD_1608FS_2; 


posit=PMD_1 608FS_1 ; 

%Convert  the  angle  to  degrees. 
cur_angle(i)=angle_calibration(h); 

probe_pos(i)=traverse_calibration(posit); 

%get  the  pressure  variable  info. 

Pressures=pressure_calibration; 
%Pressures=PMD_1 608FS_3; 
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o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%Cobra  Probe  pressures 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o 

%Right  Hole  Pressure 
P3=Pressures(1); 


%Center  Hole  Pressure 
P1=Pressures(2); 


%Left  Hole  Pressure 
P2=Pressures(3); 


%Total  Pressure  of  jet 
Poj=Pressures(4); 


%Static  Pressure  of  jet 
Psj=Pressures(5); 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%Temperature  in  pipe 

%Tp=is  entered  manually  in  the  GUI; 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%%%%%%%%%%%%CALCULATIONS%%%%%%%%%%%%%%%%%%%% 


q_probe=P1-(P2+P3)/2; 

VJet(i)=sqrt(1-(Psj/Poj)A((gam_air-1)/gam_air))*sqrt(2*Cp_air*Tp); 

Cp1_Probe(i)=(P1  -P2)/q_probe; 

Cp2_Probe(i)=(P1  -P3)/q_probe; 

Cp3_Probe(i)=(P2-P3)/q_probe; 
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At_Probe(i)=(P1  -Poj)/q_probe; 

As_Probe(i)=q_probe/(Poj-Psj); 

%Mach  Number 

Mach_No(i)=sqrt(((Poj/Psj)A((gam_air-1  )/gam_air)-1  )*(2/(gam_air-1 ))); 
end 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

%Average  the  values 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


probe_pos=(probe_pos(1 )  +probe_pos(2)+probe_pos(3))/3; 

cur_angle=(cur_angle(1)+cur_angle(2)+cur_angle(3))/3; 

Cp1_Probe=(Cp1_Probe(1)+Cp1_Probe(2)+Cp1_Probe(3))/3; 

Cp2_Probe=(Cp2_Probe(1)+Cp2_Probe(2)+Cp2_Probe(3))/3; 

Cp3_Probe=(Cp3_Probe(1)+Cp3_Probe(2)+Cp3_Probe(3))/3; 

VJet=(VJet(1)+VJet(2)+VJet(3))/3; 

At_Probe=  (At_Probe(1  )+At_Probe(2)+At_Probe(3))/3; 

As_Probe=(As_Probe(1)+As_Probe(2)+As_Probe(3))/3; 

Mach_No=(Mach_No(1)+Mach_No(2)+Mach_No(3))/3; 

Tj=Tj; 
m=1 ; 


%Output  the  data  to  the  GUI 

O/  O/  O/  O/  O/  O/  O/  O/  O/  Of  O /  Of  o /  o /  o /  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

handles  =  guihandles(gcbo);  %  generate  handles  struct 


%Set  the  pressure  values  to  the  gui  textbox 
set(handles.P_Left,  'String',  P3); 

set(handles.P_Center,  'String',  P2); 

set(handles.P_Right,  'String',  PI); 
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O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  O/  Of  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 

/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%Output  and  calculate  the  data  from  the  Three  hole  probe 

[Pitch_m,Yaw_m,Mach_m,ps_m,po_m]  =  three_hole_probe_output(m,P3,P2,P1  ,Poj); 

%Output  a  second  time  prelacing  Poj  with  P2  to  correct  for  crazy  numbers 
%near  the  casewall 

[Pitch_m_2 ,  Yaw_m_2 ,  Mach_m_2 ,  ps_m_2 ,  po_m_2] 
three_hole_probe_output(m,P3,P2,P1  ,P2); 


%flow_angle  is  the  true  angle. 
flow_angle=cur_angle+Yaw_m_2-9; 


%Create  the  data  array  to  be  written. 

header=[probe_pos,  cur_angle,  Pressures(3),Pressures(2),Pressures(1),  Pressures(4), 
Pressures(5),  Tp,  Tj,  q_probe,  VJet,  Cp1_Probe,... 

Cp2_Probe,  Cp3_Probe,  At_Probe,  As_Probe,  Mach_No, 
Pitch_m,Yaw_m,Mach_m,ps_m,po_m,  Pitch_m_2,Yaw_m_2,Mach_m_2,ps_m_2,po_m_2, 

flow_angle]; 

%Set  the  value  of  the  Mach  No  textbox  to  the  current  Machno. 
set(handles.probemach,  'String',  Mach_m_2); 


%% 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%NULL  YAW  LOOP 

%Added  to  rotate  the  probe  directly  into  the  flow  for  increased  accuracy 
%when  evaluating  the  flow  behind  the  rotor. 


%flowangle  is  the  angle  relative  to  the  probe. 
flowangle=Yaw_m_2; 
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if  abs(flowangle)<5  |  abs(flowangle)>60; 


break 

elseif  flowangle>=0; 

dur=abs(flowangle/2.9); 

ERB24_RIGHT(dur); 

pause(.5); 

current_angle; 


elseif  flowangle<=0; 
dur=abs(flowangle/3); 
ERB24_LEFT(dur); 
pause(.5); 
current_angle; 


end 

end 

o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


%% 


%Freddy.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  headers  will  be  created, 
load  freddy.txt 

%restart.txt  is  a  notepad  value  used  to  hold  on  to  data  outside  of  the 
%program  to  ensure  data  will  not  be  overwritten  if  the  run  is  restarted, 
load  restart.txt 

a=freddy+1+restart; 

%create  a  variable  that  can  be  evaluated  so  that  the  values  of  lines 
%can  be  changed. 

fred  =  ['ExcelWrite('  ""  'C:\Documents  and 

Settings\user\Desktop\Transonic\testing1  .xls . ,  0,'... 

. ['  num2str(a)  7  num2str(2) '] . ,  header)’]; 

%evaluate  the  variable  and  send  the  info  to  ExcelWrite  to  evaluate  it  and 

%send  it  to  an  EXCEL  Spreadsheet. 

eval(fred); 

%pause  to  ensure  Excel  is  closed  for  a  slow  computer 
pause(l) 

%m=1; 

%[Pitch_m,Yaw_m,Mach_m,ps_m,po_m]  =  three_hole_probe_output(m,P3,P2,P1  ,Poj) 
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%Change  the  value  of  freddy  to  next  row  number. 
freddy=freddy+1 ; 

%Save  the  row  number  to  a  text  file  so  that  it  will  not  be  cleared, 
save  freddy.txt  freddy  -ascii 


function[]  =  sampledata_run(begpoint,  endpoint,  points,  Tp,  Tj) 


%numpos  is  the  input  number  of  positions  to  take  a  sample  of. 
%range  is  the  distance  to  take  samples  on  with  zero  beginning  at  the 
%casewall 

thelength=beg  point-endpoint; 

%Create  a  linear  array 
thepos=iinspace(0,  1,  points); 


%output  an  array  of  positions  based  on  cosine 
cospos=abs(1  -cos(thepos*pi))/2; 

%create  the  array  of  final  sample  positions 
sample_pos=cospos*thelength; 

%Change  the  starting  area,  if  needed. 
flipped=fliplr(sample_pos); 


o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/  o/ 
/o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o  /o 


final_pos=flipped+endpoint; 
%plot(thepos,  final_pos) 
%break 


%for  j=1  :points 

%handles  =  guihandles(gcbo);  %  generate  handles  struct 

%set(handles.Stop_Points(j),  'String',  final_pos(j)); 

%end 
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for  i=1  :points 


Pos_lnput=final_pos(i); 

handles  =  guihandles(gcbo);  %  generate  handles  struct 

set(handles.Pos_lnput,  'String',  Posjnput); 

probe_pos_req=Pos_lnput; 

traverse(probe_pos_req); 

pause(l); 

OUTPUT_TO_EXCEL(Tp,  Tj); 
end 
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APPENDIX  O:  DIFFUSION  FACTOR  CAFCULATIONS 


The  diffusion  factor  is  a  method  of  assessing  the  blade  loading  in  an  axial 
compressor.  For  a  simple,  two-dimensional  geometry,  the  diffusion  factor  is  calculated 
using  the  equation: 

V  2  V02 

DF  :=  1 - + - 

V  i  2-g-V  i 

where  Vi  and  V2  are  the  inlet  and  exit  velocities  relative  to  the  blade,  V02  is  the 
turning  of  the  flow,  assuming  the  flow  is  axial  at  the  blade  entrance,  and  a  is  the  solidity 
of  the  blade  corresponding  to  the  flow.  It  is  believed  that  values  above  0.6  can  indicate 
blade  stall  and  values  of  0.45  may  be  a  typical  design  choice.  (Cumpsty,  1989) 

The  diffusion  factor  was  calculated  for  each  set  of  speed  data  and  results  were 
then  compared  to  CFD  data  created  with  SWIFT  Code.  A  sample  set  of  calculations  for 
100%  Speed,  Peak  Efficiency,  at  a  point  near  the  midpoint  of  the  blade  is  show  in  the 
following: 

Ds  is  the  design  speed  in  RPM. 


D  s  :=  27085 


co  is  the  angular  velocity  in  rad/sec. 


co  :=  D 


71 

60 -sec 


co  =  2.836  x  10"^Hzi 
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mdot  is  the  mass  flow  rate  in  Kg/sec. 


mdot  :=  7.687 


kg 


sec 


T0  is  the  Stagnation  Temperature,  measured  in  the  upstream  flow. 


T  0  :=  292K 

M  is  the  Mach  number,  measured  by  the  three  hole  pressure  probe  in  the 
upstream  flow.  M  has  been  area  averaged  over  the  15  data  points  taken  upstream. 


M  :=  .3344 


X  is  the  Specific  Heat  Ratio.  The  value  of  1.41  is  used  throughout  the 
calculations  for  air  and  was  used  during  the  calibration  as  well. 

X  :=  1.41 


The  Static  Temperature,  Ts  is  calculated  from  the  Mach  number  and 
Stagnation  temperature. 


T 


s 


T  s  =  285.456K. 
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Ps  is  the  Stagnation  Pressure,  measured  with  a  pitot  static  tube  upstream  of 
the  flow. 


P  s  :=  69153Pa 


Rd  is  the  Ideal  Gas  Constant  specific  for  air. 


R 


d  := 


287- 


J 

KM 


The  density,  p  is  calculated  upstream  in  the  flow,  prior  to  decrease  in  area 
caused  by  contraction  in  area  prior  to  entering  the  turbine  blades. 


P  := 


p  =  0.844 


Rd'Ts 
kg 


m 


The  diameter  of  the  inlet,  d,  is  0.2794  meters  (1 1  inches). 


d  :=  .2794 -m 
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At  is  the  total  area  of  the  inlet. 


At  =  0.061312  m2 


Knowing  the  mass  flowrate,  density,  and  area,  allows  us  to  calculate 
V upstream?  the  velocity  of  the  fluid  prior  to  reaching  the  spinner. 


mdot 

^upstream  :=  7“ 

P  'At 


m 

^upstream  =  148.534  — 


The  area  of  the  inlet  is  reduced  as  the  air  is  directed  through  the 
compressor  blades.  The  spinner  and  reduction  in  area  acts  as  a  nozzle,  and  a  new 
velocity  must  be  calculated  to  better  approximate  the  velocity  at  the  leading  edge 
of  the  blades.  The  restriction  in  area  is  treated  as  a  simple  nozzle  in  the  following 
equations. 


Entrance  •  \ 


n- (2 -.073025  -m)' 


^entrance  0.044559  m 
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The  critical  area  Astar  can  be  calculated  assuming  isentropic  flow  and  that 
there  are  no  shock  waves  in  the  duct.  (White,  2003) 


_ (AfM) 

1  +  ,5-U  -  l)-M2] 

,5-U  +  l) 


.5-(k+l)-(k-l) 


Agtar  =  0.022231m2 


Once  the  critical  area  is  know,  the  same  equation  is  then  used  to  solve  for 
the  Mach  number  at  the  entrance  to  the  compressor  blades.  The  sixth  order 
equation  was  solved  using  Matlab. 


(A  tMa) 


star 


_1  +  .5 -(A,  -  l)-Ma2] 
,5-U  +  l) 


,5-U+l)-U-l) 


Of  the  six  solutions,  only  two  were  non-imaginary.  The  first  could  not  be 
used  as  it  was  greater  than  1  and  the  second  satisfied  the  requirements  of  being 
less  than  one: 


Ma  :=  0.461935 
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The  Mach  number  was  then  designated  Mi,  the  mach  number  of  the  flow 
directly  before  the  leading  edge  of  the  compressor  blades. 

M  ]  :=  .461935 

The  density  of  the  entrance  flow  could  then  be  calculated. 


pe;=P 


1  +  ,5-U  -  l)-Mj2 


Pe 


0.937 


kg 

3 

m 


The  Static  Temperature  at  the  entrance  Tse  could  then  be  calculated,  assuming 
isentropic  conditions,  as  the  Stagnation  Temperature  T0  remains  constant  in  a  nozzle. 
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With  the  Mach  number  and  Static  Temperature,  the  Static  Velocity  Ci  could  then 
be  calculated. 

C  i  :=  M  i  d'T  se 

Ci  is  then  the  velocity  of  the  flow  as  it  enters  the  blade  row. 

m 

C  i  =  155.427  —  . 

s 

The  radial  distance,  re,  at  the  leading  edge  of  the  blade  was  used  for  calculating 
the  relative  velocity  of  the  blade  at  the  flow  entrance.  The  distance  was  calculated 
matching  the  flow  area  where  the  data  point  was  taken  downstream,  to  the  flow  area  at 
the  trailing  edge  of  the  blade.  This  flow  area  was  then  matched  to  streamlines  calculated 
by  the  CFD  Swift  code. 

r  e  :=  .1 1925m 


Entrance  Flow  Velocity  Triangle 


*Ti 


Ui  is  the  velocity  at  the  leading  edge  of  the  blade  at  the  point  that  is  measured. 
The  points  that  are  used  for  Ui  are  the  area  averaged  points  corresponding  to  the  area 


135 


averaged  points  sampled.  It  is  a  product  of  the  rotational  speed  and  radial  distance  from 
the  center. 

Uj  :=  co  -re 

m 

Ui  =  338.233  — 
s 

V i  is  the  inlet  velocity  relative  to  the  blade. 


2  2 
C1  +U1 


Vl 


372.235 


m 

s 


Vei  was  assumed  to  be  zero,  as  the  flow  was  found  to  be  axial  upstream  of  the 
compressor. 


The  flow  variables  were  then  calculated  downstream  of  the  blade  row.  The  flow 
was  sampled  at  15  positions  using  the  pressure  probe. 


M2  is  the  Mach  number  measured  of  the  flow  measured  with  the  three-hole 
pressure  probe  at  point  four  downstream  in  the  case  wall.  The  position  was  taken  near 
the  midpoint  of  the  flow. 


M2  :=  .555 
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The  Mach  number  allows  us  to  calculate  the  static  temperature,  assuming  the 
Stagnation  Temperature  has  remained  constant  through  the  compressor. 


T 


s2  := 


Ts2  =  274.657  K 


The  Static  Temperature  allows  us  to  calculate  the  Static  Velocity  of  the  flow 
downstream  of  the  compressor  rotor. 


C2  :=  M2^RdVTs2) 

C2  =  185.029  — 
z  s 

The  solidity,  a  was  matched  to  the  area  integral  of  flow  of  the  trailing  edge  of  the 
blade.  The  solidity  was  calculated  using  flow  lines  taken  at  49  points  from  a  CFD  grid 
along  the  actual  blade  dimensions  following  the  flow  line. 

The  solidity  data  was  graphed  and  is  shown  below.  Solidity  calculated  matched 
the  predicted  values  by  Sanger.  (Sanger,  1996) 
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The  solidity  was  then  interpolated  for  the  points  required  for  each  experimental 
data  point  based  on  the  ratio  of  the  area  of  flow. 

For  this  data  point  solidity  was  found  to  be: 

a  :=  1.642 


The  value  for  p,  the  yaw  angle  of  the  flow  was  provided  by  the  three-hole  probe 
from  the  values  taken  downstream  of  the  rotor.  The  angle  was  used  to  calculate  the 
velocity  relative  to  the  blade.  The  velocity  triangle  is  shown  below  in  Figure  45. 


P  :=  25.27 -deg 


138 


Exit  Flow  Velocity  Triangle 


The  variable  Cel  is  the  axial  component  of  the  static  flow  velocity  and  is 
calculated  using  the  static  velocity  and  the  angle  of  the  flow. 


C02  :=  C2  sin(p) 

Ca9  =  78.986  — 
s 

Cx2  is  the  radial  component  of  the  static  flow  velocity  and  is  also  calculated  with 
the  static  velocity  and  angle  of  the  flow  found  by  the  probe. 
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^X2  :=  C2-cos(p) 


CX2  =  167.322  — 


U2  was  the  velocity  at  the  trailing  edge  of  the  blade  at  the  area  averaged  point  of 
measurement.  Like  Ui  it  was  a  product  of  the  rotational  speed  and  the  distance  from  the 
center  of  the  measurement  point,  however,  in  this  case  the  point  was  projected  onto  the 
trailing  edge  of  the  rotor  blades.  The  radial  velocity  is  found  by  measuring  the  distance 
the  point  was  from  the  case  wall,  subtracting  it  from  the  full  radius,  and  multiplying  the 
result  by  the  angular  velocity. 


U2  :=  (0.1397  -  0.0228223)m-co 


U9  =  331.504  — 
s 

Ve;  was  the  angular  relative  velocity  and  was  the  difference  between  the  velocity 
of  the  trailing  edge  and  the  axial  component  of  flow. 


V  02  :=  U  2  “  C  02 


m 

Vfl9  =  252.518  —  1 
s 


V2  was  the  velocity  of  the  flow,  relative  to  the  blade.  With  the  angular  and  axial 
flow  components  found,  V2  was  calculated. 
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V7  =  302.923  —  1 
s 

With  all  components  of  the  equation,  the  Diffusion  Factor  can  then  be  determined 
with  the  diffusion  factor  equation: 


V  2 

DF  :=  1 - + 

VI 


V  02 
2-ct-V  j 


DF  =  0.393 . 
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